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ABSTRACT 

A  VLF  relative  navigation  system  makes  use  of  the  fact  that,  at 
any  given  point  on  the  earth,  phase  delay  of  a  received  VLF  signal  is 
highly  stable  and  predictable.   As  the  receiver  is  physically  moved, 
phase  delay  changes  linearly  with  distance  from  the  transmitting  sta- 
tion, so  that  by  keeping  track  of  the  phase  delay  of  the  received 
signal  from  several  VLF  stations,  one  may  keep  an  accurate  plot  of 
geographical  position. 

Two  problems  experienced  in  measuring  phase  delay  are  the  diurnal 
shift  and  long  path  interference.   An  investigation  is  made  into  a 
method  of  predicting  the  diurnal  phase  shift  and  the  resultant  phase 
due  to  simultaneous  reception  of  long  and  short  path  signals.   Also 
investigated  is  a  receiving  antenna  having  a  cardioid  shaped  radia- 
tion pattern  which  could  provide  discrimination  against  long  path 
signals . 
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1.   Introduction 

Certain  characteristics  of  the  Very  Low  Frequency  portion  of  the 
electromagnetic  spectrum  make  feasible  many  important  applications,  some 
relatively  new  and  others  not  so  new.   One  striking  application  is  the 
transmission  of  radio  waves  to  submerged  submarines  which  have  not  even 
an  antenna  above  the  surface.   Another  is  the  transmission  of  messages 
to  mobile  units  operating  halfway  around  the  world. 

The  application  which  inspired  this  paper,  however,  is  that  of  a 
VLF  relative  navigation  system.   Such  a  system,  using  VLF  transmitting 
stations  already  broadcasting  for  other  purposes,  has  been  tested  and 
found  feasible  for  both  ships  and  aircraft  operating  at  long  distances 
from  the  transmitting  stations.   The  system  makes  use  of  the  fact  that, 
at  any  given  point  on  the  earth,  phase  delay  of  a  received  VLF  signal  is 
highly  stable  and  predictable.   As  the  receiver  is  physically  moved, 
phase  delay  changes  linearly  with  distance  from  the  transmitting  sta- 
tion, so  that  by  keeping  track  of  the  phase  delay  of  the  received  signal 
from  several  VLF  stations,  one  may  keep  an  accurate  plot  of  geographical 
location. 

This  paper  attempts  to  resolve  one  particular  problem  associated 
with  the  relative  navigation  system  application.   This  problem  is  the 
result  of  interference  noted  from  the  reception  of  a  VLF  signal  which 
has  travelled  along  a  great  circle  path  more  than  halfway  around  the 
circumference  of  the  earth,  i.e.,  interference  of  the  so-called  long 
path  signal  with  the  short  path  signal.   Section  2  is  a  study  of  as- 
pects of  propagation  peculiar  to  the  VLF  spectrum,  while  Section  3 
discusses  in  detail  the  VLF  relative  navigation  system.   Sections  4 
and  5  provide  information  on  VLF  receivers  and  frequency  standards 


applicable  to   the  navigation  system. 

Discrimination  against   the  long  path  signal    is  desirable  to  mini- 
mize navigational   errors;    this   may  be  accomplished   easily  aboard   ship 
by  the  use  of   an  antenna   system  having   a  cardio id-type   radiation  pattern. 
VLF  antennas   are  discussed    in  Section   6,   with  emphasis  on  the  use  of   a 
loop   in  conjunction  with  a  whip  antenna  to   provide   the  desired   cardio id 
pattern. 

Sections    7   and  8   provide  a  detailed  discussion  of  the   prediction 
of   attenuation  and  phase  delay  of   VLF  waves.      These  predictions    are 
necessary   if  one   is  to   calculate   in   advance  the  effect  of   the   long  path 
signal  on  the  observed   resultant   phase  delay  at    any  given  geographical 
location.      Section  9   provides   an  example  of   just   such  a  calculation. 
Observed  data  on   long  path  signals   could  not  be  obtained   locally,    since 
there  were  no  VLF  stations   at   a  sufficient  distance  from  Monterey,    Cal- 
ifornia,   to   provide  a  measurable   long  path  signal.      It  was   thought 
that   station  GBR,   transmitting  on   16.0   Khz   from  Rugby,    England,   would 
provide  this  data;    however,    this   station  was  off   the   air  during  the 
period   available   for  experimental    investigation    (July-November   19  66). 
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2.   Characteristics  of  Very  Low  Frequency  Radio  Propagation 

The  characteristics  of  VLF  propagation  which  make  these  frequencies 
useful  for  long  distance  communication  and  highly  accurate  (to  within  a 
few  parts  in  1011)  frequency  comparison  are  a  relatively  low  path  atten- 
uation which  is  stable  with  time  and  a  phase  shift  which  is  very  constant, 
and  predictable,  over  long  periods  of  time.   VLF  also  overcomes  (to  a 
large  extent)  the  two  major  disadvantages  of  HF  transmissions,  unrelia- 
bility during  ionospheric  disturbances  and  fairly  rapid  and  deep  fading 
[28]. 

The  major  disadvantages  of  VLF  transmissions  are  the  very  narrow 
(20-150  Hz)  bandwidths  available  and  the  large  transmitting  installations 
that  are  necessary.   The  narrow  bandwidths  available  are  due  to  the  small 
spectrum  (10-30  kHz)  and,  mainly,  to  the  high  Q  of  VLF  antenna  systems. 
Since  these  antennas  are  very  short  with  respect  to  a  wavelength  their 
input  impedance  consists  of  a  large  reactance  and  a  very  small  radiation 
resistance.   To  obtain  maximum  efficiency,  therefore,  it  is  mandatory 
that  all  losses  in  the  antenna,  loading  coils,  and  ground  be  kept  to  a 
minimum.   The  resulting  low  resistance  coupled  with  high  reactance  gives 
these  systems  an  inherent  high  Q  and  consequent  narrow  bandwidth. 

The  high  powered  transmitters  are  required  to  obtain  large  areas 
of  coverage  and  because  of  the  high  level  of  atmospheric  noise  present 
at  VLF.   This  atmospheric  noise  is  due  mainly  to  lightning  discharges 
which  are  centered  in  frequency  at  about  11  kHz.   A  typical  discharge 
will  radiate  approximately  10->  joules  and  have  a  3db  bandwidth  of  about 
12  kHz  [6j. 

Theory  of  VLF  Propagation 

The  two  theories  that  have  been  predominantly  used  to  explain  the 
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characteristics  of  VLF  propagation  are  the  ray  theory  and  the  waveguide 
mode  theory.   It  has  been  stated  £28}  that  the  mode  and  ray  theories  are 
related  and,  if  all  modes  and  all  rays  were  taken  into  account  with  ap- 
propriate corrections  for  losses,  both  methods  would  yield  the  same  re- 
sults.  However,  at  long  distances  the  mode  theory  is  simpler  to  use 
because  only  one  mode  need  be  considered  but  many  individual  rays  must 
be  taken  into  consideration.   At  short  distances  many  modes  are  present 
making  the  mode  theory  calculations  complicated,  but  only  a  few  rays  need 
be  considered  since  the  higher  order  rays  are  incident  on  the  ionosphere 
at  much  smaller  angles  where  the  reflection  coefficient  is  small.   In 
this  paper  the  only  distances  which  will  be  considered  are  those  for 
which  the  mode  theory  calculations  may  be  assumed  valid  with  only  one 
mode  present  (i.e.,  greater  than  4000  km  from  the  transmitter). 

From  observed  measurements  of  the  electric  field  within  the  earth- 
ionosphere  waveguide  it  has  been  found  that  the  model  shown  in  Figure 
2-1  can  be  used  to  approximate  the  actual  physical  case.   The  mode  of 
principal  interest  at  long  distances  is  the  TM, ,  shown  in  the  left  side 
of  Figure  2-1,  the  TM~  and  higher  order  modes  being  small  enough  in 
amplitude  beyond  4  Mm  to  be  neglected.   From  Watt  and  Croghan  Cl83>  for 
distances  greater  than  1  Mm  the  vertical  electric  field  can  be  represen- 
ted as  the  sum  of  a  very  few  waveguide  modes: 

n  =  1  to  3 
Ez  =  J2  EZjn      d>l  Mm 

n  =  1 

Here  E    is  the  electric  field  strength  for  each  particular  mode   (in 

db  relative  to  1  v/m) 


Ez  n  [db,  v/m]  =  10  Log  Pr  +  K^  -  10  Log  f 
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Figure  2-1. 
Electrical  field  structure  for  waveguide  modes. 

(Wait  CllD) 


13 


-  10   Log  hijt   +    10   LogAnjt 

-  10   Log    (a  sin  d/a)    -   10  Log  h-    „ 
+    10   LogAn>r   +    20   Log  Gh>n 

"   Ld,n   -°Cd/1°6 
where 

P     =  power   radiated   into  the  half   space   above  the  earth    (watts). 

Kn  =   104.3  db   for  modes    1,2,    a   constant  which  relates   field 

strengths   at  the  surface  to   power   radiated. 

f   =  frequency   (Hz) 

h-  +.  =  effective  ionospheric  height  at  the  transmitter  (m) . 

h.    =  effective  ionospheric  height  at  the  receiver  (m). 
l  ,r 

JL     .    =  mode  excitation  modification  factor,  defined  as  the  ratio 
of  power  launched  into  the  concentric  spherical  shell  guide  relative  to 
that  for  a  flat  guide  with  perfectly  conducting  earth, 
a  =  earth's  radius  (»6.4  Mm) 
d  =  path  length  (m) 

(a  sin  d/a)  accounts  for  energy  spreading  in  the  spherical 
guide  structure. 

An  r  =  mode  field  modification  factor  for  a  curved  earth  at  the 
ionospheric  height  at  the  receiver. 

20  Log  G^   =  height  gain  factor  (db)  relative  to  the  field  at 
the  surface. 

L,   =  loss  due  to  discontinuities  along  the  path  (db)  . 

O1^   -   effective  attenuation  rate  in  db/Mm  (different  for  each 
mode) . 

These  factors  are  discussed  at  some  length  by  Watt  and  Croghan  and  will 
not  be  commented  on  further  in  this  paper,  with  the  exception  of  the 
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attenuation  factor. 

For  paths  with  uniform  boundary  surfaces,  the  attenuation  factor  is 
given  by 

OC  =  M(0)Ki  +  CXg 
where 

M( 0 )    =   earth's  magnetic   field   factor 

_        attenuation  with  magnetic   field 
attenuation  without  magnetic   field 

<P  =  direction  of   propagation  with  respect   to   the  earth's 
magnetic   field 

CX-    =   attenuation  constant  of   the   ionosphere 
0<g   =   attenuation  constant   of  the  ground. 
CX£   and  CXg   are  related  to   the  reflection  coefficients  of  the   earth  and 
the   ionosphere  which  are   in  turn   related  to   the  conductivities  of  these 
two  media. 

The   attenuation  curves    shown   in  Figures    2-2   and   2-3  are  derived 
for  M(0)    =    1,    i.e.    independent  of   the  earth's  magnetic   field.      Here 
X    is   the   conductivity  gradient  of   the   ionosphere.      From  these  curves   the 
difference    in   attenuation  of    the   first   two  modes    is   apparent.      At    15  kHz 
and   for   a  4  Mm  path  the   smallest  difference   in   signal   strength  between 
the  two  modes    is  about    13  db,    justifying  the  earlier  assumption  that 
only  the   first  mode  need   be  considered   at   these  distances.      Figures 
2-4,    2-5,    and  2-6  are   shown  for  comparison  and  while  there   are   some  dif- 
ferences   in  the  attenuation  curves   for  n  =    1,    the  n  =   2   curves   show  the 
same  larger   attenuation  rates    (with  respect   to   the   n  =   1   curves).      Fig- 
ure 2-7,    while  not   a  set  of   calculated  curves,    shows   very  vividly  the 
effect  of   the  earth's  magnetic   field  on  the   attenuation  factor. 
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The  phase  velocity  of  a  certain  mode  for  a  fixed  set  of  conditions 
can  be  determined  from  a  solution  of  the  mode  equation.   This  is  done 
by  Wait  Cn!3  an^  Wait  &  Spies  [_23]    resulting  in  the  curves  shown  in 
Figures  2-8  through  2-11.   The  effect  on  phase  velocity  of  several  fac- 
tors can  be  seen  here.   Assuming  a  constant  frequency,  it  can  be  seen 
that  the  height  of  the  ionosphere  has  a  large  effect  on  the  phase  velocity, 
lowering  it  with  increasing  height.   This  is  the  cause  of  the  familiar 
diurnal  shift  of  phase  delay  of  a  received  VLF  signal.   At  night  the  D 
layer  of  the  ionosphere  (the  lowest  layer)  dissipates  and  the  effective 
height  of  the  ionosphere  is  increased  by  about  20  km.   The  resultant  de- 
crease in  phase  velocity  causes  an  increase  in  the  phase  delay  of  the 
received  signal  as  compared  with  a  local  standard.   After  sunrise  the 
sun's  rays  ionize  the  lower  portions  of  the  ionosphere,  the  D  layer  re- 
turns, and  the  process  repeats  itself. 

Comparing  Figures  2-8  and  2-9  will  show  that  an  imperfectly  re- 
flecting ionosphere  (e.g.  ,  an  exponential  variation  of  ion  density)  will 
increase  the  phase  velocity  slightly.   The  effect  of  finite  ground  con- 
ductivity is  seen  in  Figure  2-10,  indicating  that  the  phase  velocity 
will  be  slightly  lower  over  land  (as  compared  to  sea).   From  Figure  2-11 
it  appears  that  the  effect  of  the  earth's  magnetic  field  on  phase  vel- 
ocity is  small  in  comparison  to  the  other  factors.   (The  effect  is  much 
more  pronounced  for  the  n  =  2  mode;  however,  that  mode  is  being  consid- 
ered negligible  in  amplitude  here).   Wait  and  Spies  have  many  more  curves 
of  phase  velocity,  attenuation  factors,  excitation  factors,  etc.,  which 
would  be  of  interest  to  anyone  planning  a  more  detailed  study  of  VLF 
propagation. 
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Experimentally  Obtained  Data 

Figures  2-12  through  2-14  along  with  the  experimental  curves  in 
Figures  2-2,  2-3  and  2-7,  are  examples  of  the  vast  amount  of  experi- 
mental data  taken  on  the  attenuation  rates  and  phase  velocities  of  VLF 
radio  waves.   The  agreement  of  this  data  with  calculated  curves  in  cer- 
tain cases  (e.g.,  Figures  2-13  and  2-14)  is  remarkable  and  provides  a 
basis  for  using  the  calculated  curves  to  predict  phase  velocities  and 
attenuation  rates. 

Figure  2-12  shows  the  many  factors  that  must  be  taken  into  account 
in  predicting  attenuation  rates:  time  of  day,  direction  of  propagation, 
frequency,  and  ground  conductivity.   In  the  10-30  kHz  frequency  range, 
the  range  of  attenuation  rates  is  about  4  db/Mm  (minimum  to  maximum) 
indicating  that  for  long  distances  a  considerable  difference  in  signal 
strength  could  be  expected  depending  on  the  path  traveled. 

In  this  paper  the  experimental  data  will  be  used  wherever  possible 
and  the  calculated  curves  will  be  used  as  a  check  and  to  fill  in  the 
gaps  not  covered  by  experimental  evidence  (e.g.,  attenuation  of  the 
n  =  2  mode) . 

Summary 

The  one  characteristic  of  VLF  radio  transmissions  of  interest  in 
this  paper  is  the  good  phase  stability  of  the  signals  over  long  periods 
of  time  and  the  predictability  of  this  phase.   Because  of  this  phase 
stability  frequency  comparisons  over  long  distances  are  possible  with  an 
accuracy  on  the  order  of  one  part  in  10   .   The  amplitude  and  phase  of 
these  signals  can  be  predicted  fairly  accurately  at  distances  greater 
than  one  megameter  from  the  transmitter  by  using  the  mode  theory  of 
propagation.   The  actual  solution  of  the  mode  theory  equations  is  a 
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rather  complex  problem;  however,  the  solutions  of  Wait  and  others  have 
been  published  and  show  a  very  good  agreement  with  experimental  evidence 
in  some  cases. 
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Figure  2-2 

Attenuation  rates  over  predominantly  sea  water 
paths  for  n=l  (experimental  and  theoretical). 

(Watt  &  Croghan  [ 18] ) 
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Figure  2-3 

Attenuation  rates  over  predominantly  sea  water 
paths  for.n=2  (experimental  and  theoretical). 

(Watt  &.  Croghan  £1*0) 
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Figure  2-4 

Attenuation  rate  of  the  first 
node  for  a  perfectly  conduct- 
ing earth  and  an  imperfectly 
reflecting  ionosphere, 

(Wait  Cll]) 
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Figure  2-5 
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Attenuation  rate  of  the  second 
mode  for  a  perfectly  conducting 
earth  and  an  Imperfectly  re- 
flecting ionosphere. 

(Wait  [11]) 


20 


i        I        I        i        i        i        I L 

12  '  16  20~'  '   "    24 

FREQUENCY,    kc/S  ' 


Figure  2-6 

attenuation  rate  of  the  first  mode  for  an  imperfectly 
conducting  earth  and  an  imperfectly  reflecting  ionosphere 

(Wait  till) 


Figure  2-7 

Attenuation  rates  for 
predominantly  sea  water 
paths  for  moderate  and 
low  sunspot  numbers  (ex- 
perimental) - 

(Watt  &  Croghan  Cl8]) 
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Phase  velocity  of  the  first 
or  dominant  mode  for  a  per- 
fectly conducting  earth  and 
a  perfectly  reflecting  iono 
sphere * 
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Figure  2-9 

Phase  velocity  of  the  first  mode 
for  a  perfectly  conducting  earth 
and  an  imperfectly  reflecting 
ionosphere - 

(Wait  CUD 
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Figure  2-10 

Phase  velocity  of  the  first  mode  for  an  imperfectly 
conducting  earth  an  axr  imperfectly  reflecting  iono- 
sphere - 

(Wait  £11]) 
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Figure  2-11 

Phase  velocity  as  a  function  of  frequency  showing  effect 
of  reflecting  height  and  magnetic  field  - 

' (Wait  &  Spies  [23]) 
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Figure"  2-12 


Attenuation  rates  deduced 
from  various  experimental 

data 
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Figure  2rl3 


Experimental  phase  velo- 
cities v/ith  pertinent 
theoretical  data ~ 
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Figure  2-14 


Experimental  phase  velo- 
cities with  pertinent 
theoretical  data- 

(Wait  &  Spies  {££]) 
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3.   A  Proposed  Very  Low  Frequency  Relative  Navigation  System 

The  proposed  very  low  frequency  (VLF)  relative  navigation  system 
has  been  outlined  by  Stanbrough  and  Keilly  [l6l,  and  presented  in  de- 
tail by  Lake  [25].   The  system  has  many  advantages,  providing  (1)  rela- 
tive simplicity;  (2)  a  long  range,  accurate  system  with  worldwide 
coverage;  (3)  the  capability  of  using  VLF  phase-stabilized  stations 
already  on  the  air,  with  no  requirement  for  additional  stations;  and 
(4)  reliable  24-hour  coverage.   The  only  required  equipment  for  ship- 
board use  of  this  system  are  one  or  more  phase  tracking  receivers  and, 
preferably,  a  directional  antenna  system  such  as  a  loop  used  in  con- 
junction with  a  sense  antenna  (whip).   The  latter  is  discussed  elsewhere 
in  this  paper.   One  recorder  channel  to  provide  a  continuous  record  of 
phase  for  each  tracking  receiver  used  would  also  be  helpful. 

The  extreme  phase  stability  of  VLF  radio  waves  propagated  over  great 
distances,  plus  or  minus  one  microsecond  under  normal  conditions,  is 
the  foundation  of  this  navigation  system.   An  oscillator  at  the  distant 
VLF  transmitter  provides  the  station  with  phase-stabilized  emissions, 
while  a  stable  oscillator  aboard  ship  at  the  same  frequency  provides 
the  local  reference.   This  local  oscillator  is  used  in  conjunction  with 
a  phase  tracking  receiver  to  provide  an  indication  of  relative  phase  of 
the  VLF  wave  at  the  receiver.   A  high-quality  quartz  crystal  oscillator 
would  be  suitable  for  use  aboard  ship  and  a  typical  unit  in  use  at  the 
U.S.  Naval  Postgraduate  School  provides  a  frequency  drift  (due  to  crystal 
aging)  of  about  2.5  parts  in  10   per  day.   This  drift  would  be  deter- 
mined for  the  particular  oscillator  in  use  aboard  ship  and  its  frequency 
would  then  be  corrected  on  a  regular  basis,  perhaps  daily,  to  keep  it 
within  some  specified  frequency  offset.   At  the  time  of  taking  a  fix, 
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the  phase  variation  due  to  the  frequency  offset  of  the  oscillator  would 
then  be  known  and  a  compensation  could  be  made. 

Q 

Assuming  a  nominal  value  of  3  x  10  m/sec  for  the  phase  velocity 
of  a  VLF  wave,  the  number  of  microseconds  required  for  the  wave  to 
travel  one  nautical  mile  is  found  to  be  6.18.   In  view  of  this,  it  is 
clear  that  for  each  increase  in  phase  delay  of  6.18  usee,    the  ship  has 
moved  one  nautical  mile  farther  from  the  transmitter. 

The  VLF  relative  navigation  system  makes  use  of  the  foregoing;  by 
keeping  track  of  the  phase  of  two  or  more  stations  as  the  ship  moves 
away  from  an  initial  known  point,  an  indication  of  the  ship's  geograph- 
ical position  is  available  at  all  times.   By  programming  a  computer  to 
make  the  calculations  and  necessary  corrections,  one  could  even  obtain 
a  continuous  track  of  the  ship's  movement.   It  is,  however,  necessary 
to  use  VLF  stations  far  enough  away  that  mode  interference  is  not  ex- 
perienced, and  only  the  first  waveguide  mode  is  involved  (greater  than 
4000  km). 

One  problem  complicating  the  navigation  system  is  that  of  the 
diurnal  phase  variation  experienced  with  VLF  waves.   Fortunately,  how- 
ever, this  variation  is  predictable  and  the  shipboard  navigator  can  be 
provided  with  tables  to  correct  for  the  phenomenon.   Figure  3-1  illus- 
trates the  diurnal  variation,  and  is  an  actual  plot  of  the  phase  of  the 
13.6  kHz  signal  from  the  Omega  Navigation  System  transmitter  located  at 
Haiku,  Hawaii,  referred  to  a  rubidium  vapor  frequency  standard  at  the 
U.S.  Naval  Postgraduate  School. 

A  first  approximation  to  a  solution  of  the  problem  of  diurnal 
phase  change  may  be  made  by  assuming  a  trapezoidal  pattern,  or  trapezium, 
as  shown  in  Figure  3-2.   A  comparison  with  Figure  3-1  shows  the  similar- 
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Figure  3-1.   Relative  Phase  of  Haiku  Omega  Station, 
13.6  kHz,  as  Received  at  Monterey 
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Figure  3-2.   Trapezium  Approximation  of  Diurnal  Shift, 
ity  of  the  two  patterns.   Time  intervals  AB  and  EF  have  the  same  values 
of  relative  phase  and  represent  the  period  of  time  during  which  the 
entire  propagation  path  is  in  sunlight.   Assuming  the  transmitter  loca- 
tion west  of  the  receiver,  point  B  is  the  time  of  local  (receiver)  sunset, 
point  C  indicates  sunset  at  the  transmitter,  point  D  local  sunrise,  and 
point  E  sunrise  at  the  transmitter.   It  has  been  found  that  correction 
of  sunrise  and  sunset  times  to  those  viewed  from  the  effective  height  of 
the  lower  boundary  of  the  D  layer  of  the  ionosphere  (about  70  km  during 
the  day  and  90  km  at  night)  provides  better  correlation  with  the  observed 
discontinuities  in  the  phase  plot. 

The  diurnal  shift  (At  in  Figure  3-2)  varies  with  distance  from  the 
transmitter,  since  it  is  caused  by  the  fact  that  the  phase  velocity  of 
a  VLF  wave  has  one  value  during  the  day  and  a  lower  value  at  night.   As 
discussed  elsewhere  in  this  paper,  these  phase  velocities  may  be  pre- 
dicted for  given  paths,  allowing  the  calculation  of  the  diurnal  shift 
at  any  given  geographical  position. 

Other,  more  precise,  models  of  the  phase  plot  may  be  determined 
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(as  shown  later  in  this  paper).   Whichever  model  is  used,  the  procedure 
to  be  used  in  navigation  is  to  predict  the  deviation  of  phase  from  the 
reference  (all-daylight  path)  value,  for  the  particular  time  of  the 
desired  fix,  and  subtract  from  indicated  phase.   As  mentioned  earlier, 
these  deviations  could  be  prepared  and  published  in  advance. 

An  electronic  navigation  system  such  as  that  being  considered  is 
known  as  a  circular  grid  system,  as  opposed  to  a  hyperbolic  system  such 
as  LORAN,  which  makes  use  of  phase  difference  between  a  master  and  a 
slave  station.   In  the  VLF  relative  system,  navigational  charts  may  be 
overprinted  with  circles  representing  lines  of  constant  time  delay  drawn 
about  each  VLF  station  as  a  center.   These  would  be  constructed  using 
the  6. 18 yusec/mile  discussed  earlier.   (It  should  be  noted  that  use  of 
these  charts  introduces  a  small  error  since  the  phase  velocity  of  a  VLF 
wave  is  not  exactly  equal  to  the  velocity  of  light;  hence,  actual  phase 
delay  differs  slightly  from  6. 18  //.sec/mile. ) 

In  order  to  use  such  an  overprinted  chart,  one  would  first  need  to 
plot  the  ship's  position  from  other  information,  which  would  then  be 
used  as  a  reference  for  future  fixes  obtained  by  the  VLF  relative  sys- 
tem.  The  phase  indicated  by  the  tracking  receivers  for  each  VLF  station 
being  tracked  would  also  be  noted  at  the  time  the  reference  fix  was  ob- 
tained.  At  the  time  the  VLF  fix  is  desired,  the  phase  readings  on  each 
VLF  station  would  again  be  noted.   For  those  stations  whose  propagation 
paths  are  all  in  daylight,  the  phase  difference  could  then  be  transferred 
directly  to  the  chart.   This  would  be  in  the  form  of  a  time  difference 
for  each  station  which  would  indicate  the  distance  travelled  away  from 
or  toward  the  station.   A  short  segment  of  a  circle  for  each  station 
would  be  interpolated  on  the  chart  between  those  segments  already  over- 
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printed.      The   intersection  of  these   circular   segments  would  then  be  the 
ship's   position. 

For  those   stations  whose  propagation  paths   are   partially  or  totally 
in  darkness,    the  correction   for  diurnal   shift   corresponding  to   the   ship's 
approximate  position  and  time  of   day  would  have   to   be  applied.      There 
would  be   a  correction   for  the  time  of  the   reference  position  and  one  for 
the  time  of   the   new  position.      The  corrections   would   be   subtracted   in 
both  cases,    since   phase  delay   is   greater   at   night.      After  the  diurnal 
correction  has   been  applied,   the  phase   information  is   used   as  described 
in  the   last   paragraph. 

It  is  important  to  note  that,  since  stations  being  used  for  navi- 
gation will  be  at  least  4000  km  distant,  the  circles  of  constant  phase 
will  appear  as  straight  lines  in  small  areas.  In  fact,  if  overprinted 
charts  are  not  available,  straight-line  approximations  of  the  circular 
segments  may  be  used  if  the  true  bearings  are  known  from  the  ship  to 
the  transmitting  stations.  These  could  be  obtained  from  a  previously- 
prepared  table.  The  straight-line  approximations  will  be  adequate  if 
the  ship's  position  is  known  within  about  50  miles  of  its  actual  posi- 
tion [25]  . 

A  very  complete   example  of   a   navigational   fix  obtained  by  this 
method,    and  using  the  trapezoidal   approximation  of   diurnal   phase   shift, 
is   given  by  Lake    \J2.5~} . 

There   are   certain  problems   associated  with  the  VLF  relative  navi- 
gation method,    some  of  which  are   (1)    cycle-slipping,    (2)   interference 
from  the   long  path    (around   the  world)    signal,    (3)    inexactness  of  the 
trapezoidal   model   of  diurnal   phase  change,    (4)   transmitter  outages, 
(5)    frequency  offset  of  the  transmitter  master  oscillator,    and   (6) 
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sudden  ionospheric  disturbances  (SID's). 

Cycle-slipping  is  a  phenomenon  in  which  the  tracking  receiver  indi- 
cates and  tracks  one  or  more  full  cycles  ahead  or  behind  the  true  rela- 
tive phase  of  the  VLF  wave.   This  often  occurs  near  sunrise  or  sunset 
when  relative  phase  is  changing  rapidly.   Frequently  the  diurnal  phase 
shift  will  be  more  than  one  cycle,  so  that  cycle-slipping  is  especially 
likely  to  occur.   Cycle-slipping  may  also  be  caused  by  interference  from 
waveguide  modes  higher  than  the  first  order,  and  by  interference  from  the 
long  path  signal  [27]  . 

The  long  path  signal  may  be  thought  of  as  a  phasor,  having  some 
magnitude  and  phase  relative  to  the  short  path  signal.   Together,  the 
two  signals  form  a  resultant  which  constitutes  the  composite  signal  re- 
ceived by  an  omnidirectional  antenna.   Figure  3-3  shows  the  phasor  addi- 
tion of  the  two  signals,  with  E^  =  |eJ  /0°  representing  the  short  path 
signal  (taken  with  phase  =  0°  for  convenience),  E2  =  IE2I  /9°  the  long 
path  signal,  and  ER  =  [Ep[  /if   the  resultant.   The  notation  E  is  used  to 
indicate  a  phasor  quantity.   Since  the  long  path  signal  may  travel 


Figure  3-3.   Phasor  Diagram 
hundreds  or  thousands  of  miles  farther  than  the  short  path  signal,  the 
angle  9  may  take  on  any  value  between  0°  and  360°.   It  is  clear,  then, 
that  for  any  value  of  lE^  which  is  an  appreciable  fraction  of  |E]J  ,  the 
resultant  will  vary  considerably  from  the  short  path  signal.   Since  such 
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a  condition  will  occur  under  favorable  circumstances,  the  need  for  a  di- 
rectional antenna  which  will  discriminate  against  the  long  path  signal 
is  evident.   The  long  path  signal  is  particularly  likely  to  be  appreciable 
near  the  antipode  of  the  transmitter,  as  well  as  at  other  locations  where 
the  short  path  signal  for  any  reason  undergoes  high  attenuation  concur- 
rent with  low  attenuation  of  the  long  path  signal. 

The  trapezoidal  model  of  diurnal  phase  change,  discussed  earlier  in 
this  section,  is  a  good  approximation  in  many  cases;  however,  there  are 
also  frequent  instances  in  which  it  introduces  excessive  error.   A  more 
refined  model  of  diurnal  phase  change,  useful  in  all  cases,  is  presented 
in  a  later  section. 

Transmitter  outages  may  occur  on  a  scheduled  basis  for  maintenance, 
in  addition  to  those  unscheduled  outages  due  to  equipment  failure  which 
invariably  occur.   During  the  interval  of  a  transmitter  outage,  the 
indicated  phase  of  a  tracking  receiver  remains  constant.   Since  the 
transmitter  master  oscillator  normally  continues  running  during  the 
outage,  the  receiver  will  resume  tracking  on  the  correct  cycle  when  VLF 
transmissions  resume,  if  the  receiver  has  remained  geographically  sta- 
tionary.  If,  however,  the  ship  has  moved  a  distance  on  the  order  of 
some  multiple  of  a  wavelength,  the  receiver  will  track  the  incorrect 
cycle  (the  same  one  as  if  it  had  remained  stationary) .   The  problem  of 
determining  the  correct  cycle  is  known  as  lane  identification,  and  is 
minimized  in  other  systems  by  utilizing  lower  frequency  (longer  wave- 
length) signals.   Additional  information,  such  as  the  ship's  dead 
reckoning  track  or  positional  information  from  other  means,  is  necessary 
to  identify  the  correct  lane,  or  cycle. 

Frequency  offset  of  the  transmitter  master  oscillator  of  a  VLF 
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station  is  occasionally  made,  for  example  to  bring  the  transmitted  fre- 
quency into  coincidence  with  the  United  States  Frequency  Standard.   These 
frequency  corrections  would  occur  infrequently  and  advance  notice  thereof 
could  be  given  to  users  of  a  VLF  relative  navigation  system.   With  such 
advance  notice,  the  navigator  would  have  no  difficulty  applying  a  cor- 
rection for  the  offset. 

Sudden  ionospheric  disturbances  are  caused  by  three  phenomena:   solar 
flares,  magnetic  storms,  and  high  altitude  nuclear  bursts.   While  these 
disturbances  occur  infrequently,  they  cause  phase  disturbances  on  VLF 
circuits  over  wide  areas  of  the  earth,  causing  errors  in  VLF  navigation 
in  the  case  of  severe  disturbances.   Chilton,  Steele,  and  Crombie  [193 
have  pointed  out  that  solar  flares  cause  a  depression  in  the  effective 
height  of  the  ionosphere,  which  in  turn  causes  increased  phase  velocity, 
resulting  in  a  lower  value  of  relative  phase  at  the  receiving  point.   A 
typical  flare  has  an  onset  time  of  one  to  30  minutes  and  requires  from 
one  to  two  or  three  hours  to  return  to  normal.   A  typical  phase  change 
caused  by  the  flare  might  be  10  microseconds,  but  varies  with  path 
length,  frequency,  and  amount  of  ionosphere  height  change. 
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4.   VLF  Receivers. 

The  receiver  used  to  obtain  experimental  data  on  phase  and  amplitude 
of  various  VLF  stations  for  this  paper  was  a  Tracor  Model  599G  VLF  phase 
tracking  receiver.   A  simplified  block  diagram  showing  the  essential  ele- 
ments of  such  a  tracking  receiver  is  seen  in  Figure  4-1.   The  output  of 
the  block  marked  VLF  Receiver  is  at  an  intermediate  frequency  of  1  kHz 
with  a  certain  phase  0  VLF'  anc*  *"s  one  ^nPut  °f  tne  Phase  comparator. 
The  other  input  of  the  phase  comparator,  a  1  kHz  signal  of  phase  0  »> 
is  derived  from  the  external  frequency  standard.   The  output  of  the  phase 
comparator  provides  the  electronic  phase  shifter  with  an  error  signal 
A  0  ,   which  tells  the  phase  shifter  how  much  to  shift  the  phase  of  the 
external  standard  to  maintain  phase  tracking.   The  output  of  the  phase 
shifter  then  goes  to  the  frequency  synthesizer  which  provides  the  ref- 
erence signal  (?  A.   The  condition  for  phase  tracking  is  A  0  =  0^  „  - 

0TA    =    0. 

The  phase  tracking   servo   and   long   integration  time  of   the  receiver 
provide   the  capability  of   tracking  with  an   input    signal-to-noise  ratio 
of   -50  db.      The  rated   sensitivity   is  0.01  microvolts  of   signal   at  the 
antenna  terminals  to   provide   phase   tracking.      Using   a   loop   antenna  made 
by  Tracor  especially  for  this   receiver    (Model   599-600),   this   corresponds 
to   a  signal    strength  of  0.3  microvolts/meter.      With  a   150   sec   servo   time 
constant    (front   panel    switch  provides   time  constants  of   5,    15,    50,    and 
150   sec)    the   equivalent   noise  bandwidth   is  0.002  Hz.      A  front-panel 
switch  on  the   receiver  allows   selection  of  either  a  broadband  RF  filter 
or  any  one  of   five  narrowband  ones.      The  narrowband   filters    provide  an 
RF  bandwidth  of   500  Hz,   while  the   IF  bandwidth   is    50  Hz.      The  receiver 
is   tunable   from  8.0   to    31.9  kHz    in  0.1   kHz   steps. 
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The  receiver  has  available  several  outputs  that  proved  useful  in 
experimental  work  for  this  paper:  (a)  100  usee  and  10  yUsec  full  scale 
recorder  outputs,  (b)  a  front  panel  bi-directional  digital  counter  for 
cumulative  phase  change,  (c)  a  100  kHz  coherent  output,  and  (d)  an  AGC- 
derived  relative  carrier  level  recorder  output.   The  100  usee  and  10 
jjsec   full  scale  recorder  outputs  use  the  receiver's  internal  phase  com- 
parator to  give  a  dc  voltage  level  proportional  to  the  phase  difference 
between  the  external  frequency  standard  and  the  received  VLF  signal. 
(Note:  the  phase  comparator  of  the  receiver  does  not  actually  operate 
at  1  kHz  as  indicated  in  the  discussion  of  Figure  4-1.   The  receiver 
contains  two  linear  phase  comparators;  one  operates  at  100  kHz  to  pro- 
vide the  10  ^tsec  full  scale  recorder  output,  while  the  other  operates  at 
10  kHz  for  the  100  yLtsec  full  scale  recorder  output.)   These  recorder 
outputs  were  convenient  and  proved  quite  satisfactory  for  continuous 
monitoring  of  the  phase  of  any  phase-stable  VLF  station.   Moseley  Model 
51P0B  strip  chart  recorders  were  used  for  monitoring  these  and  the 
other  receiver  outputs,  as  desired. 

The  electronic  phase  shifter  of  the  receiver  actually  operates  at 
10  MHz.   One  cycle  of  phase  shift  then  corresponds  to  0.1  usee  drift 
which  feeds  the  front  panel  bi-directional  phase  counter.   This  digital 
counter  provides  a  cumulative  indication,  so  that  by  taking  readings  at 
two  different  times,  one  may  determine  by  subtraction  the  exact  phase 
drift  that  has  occurred  during  the  elapsed  time. 

The  100  kHz  coherent  output  of  the  receiver  provides  a  100  kHz 
signal  from  the  frequency  synthesizer  section  shown  in  Figure  4-1.   This 
signal  thus  has  a  phase  corresponding  to  that  of  the  received  VLF  signal 
This  output  proved  useful  as  an  input  to  an  external  phase  comparator, 


37 


allowing  a   check  on  the   phase  of   any   local    frequency   standard   (the  other 
input  to   the   comparator)   with  respect   to   the  received  VLF  signal.      Both 
the  rubidium  vapor   frequency   standard   and  the  quartz  crystal   oscillator 
in  use   at   the  Naval    Postgraduate   School  have   100  kHz  outputs   suitable 
for  this    (see   Frequency   Standards   section). 

The  relative  carrier   level  output  of  the  tracking  receiver  provides 
a  dc    level   proportional   to  the  AGC  voltage  developed   in  the   receiver   and 
is  thus   a   relative   indication  of   received  signal   strength  of   the  VLF 
station.      The  AGC  has    a  range  of   approximately  kO  db,    and   peaks   at   about 
40  db   above  the  minimum  signal   necessary  for  tracking. 

One  problem  occasionally  encountered   in  experimental   work  was   the 
very   long   integration  time  of   the  receiver.      For  example,    a  cardioid 
antenna  system  must   be  adjusted   for   proper   phase  and   amplitude  of   the 
whip  antenna  when  changing  to   a  new  frequency.      To  get   a  good   indication 
of   phase   and  amplitude,    it  was   usually  necessary   to  wait  five   to    10  min- 
utes  after    each  adjustment  of   the  whip.      Since  several   adjustments   are 
normally  required,   the  process   can  be   quite  time-consuming. 

Adjustments   are  particularly  time-consuming   if  one   is  receiving 
one  of  the  Omega  stations.      Each  of   these  stations   transmits  only  about 
one   second  out   of   10,    so   the   long   integration  time   problem   is   compounded. 
At  the  time  of   the  experimental  work   for  this   paper,    there  were  four 
Omega  stations   in  operation    (there  are  planned   ultimately  to   be  eight). 
Each  station  operates   on  a  time-shared  basis  on  two   frequencies,    10.2 
and   13.6  kHz.      They  transmit   an  unmodulated  waveform   in  a  prescribed 
sequence   for   a  prescribed   length  of   time    (e.g.,   0.9   sec,    1.1   sec,    etc.) 
to  facilitate  identification.      Still    in  the  planning  stage   for  the  sys- 
tem are  the   addition  of    11.33  kHz   by   all   stations  on  the  time-shared 
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basis,   as   well   as   transmissions  of   12.00  and   12.75   kHz    (not   on  a  time- 
shared  basis)   by   some  of   the  stations.      In  order   to  obtain  phase  and 
amplitude   information  on  an   individual    station,    it    is   necessary  to 
clamp  the  AGC   and  blank  the  receiver  during  the  time  the  desired   station 
is    not  transmitting.      In  the   laboratory  at   the  Naval   Postgraduate   School, 
this  was   accomplished  by  a  Tracor   Omega  Gating  Unit.      This   unit  may  be 
used  with  four  tracking   receivers   simultaneously,    and   a  modification  for 
the  receiver   is   provided  which  allows   switching  between   normal   receiver 
operation  and  Omega  operation.      The  unit   proved  highly   satisfactory   and 
much  valuable  data  was   obtained  using  Omega  stations. 

It  was   useful   in  the   experimental   work  to  have  an  auxiliary  VLF 
receiver   available  which  could  conveniently  be  tuned   from  station  to 
station  for  quick  audio    identification.      The  receiver  used  was   a  model 
AN/SRR-11,    and  a  similar  type  would   no   doubt  be  useful   for  auxiliary- 
purposes    in  the   shipboard   VLF  relative  navigation  system. 
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5.      Frequency  Standards 

A  good   frequency   standard  must   be  used   at   both  the  phase-stabilized 

VLF  station  and    the   ship   using   the  VLF   relative   navigation   system;    hence 

it    is   a  vital   part   of  the   system.      Its   use   in  the  system  aboard   ship   is 

« 
shown   in  Figure   4-1. 

There  are   two   basic   types   of   precise   frequency   standards:    the  quartz 
crystal  oscillator   and  the  atomic   frequency  standard.      The   former   is 
characterized  by  a   predictable    frequency  drift   as  the  crystal   ages, 
while  the   latter  has   either   an  extremely  small   drift  or  none  at   all. 

The   past   twenty  years   have  seen  great    achievements    in  the  field 
of   precise  frequency  control;    today,   time  difference  or   frequency  may 
be  measured  more  precisely  than  any  other  physical   quantity.      Using  an 
atomic   hydrogen  maser,    a  frequency   stability    (RMS   deviation  from  the 
mean)   of   2   parts    in  10      ,   over  a  one-day  averaging   period,   may  now  be 
obtained    [30]  .      One  problem  which  has    arisen   in  the  growth  of   technol- 
ogy of   frequency   standards    is   a   lack  of    standardization  of  the  defini- 
tions  of   certain  terms.      Two   of  these  are   short-term  stability  and 
long-term  stability.     The  former  refers   to   average   frequency  deviation 
measured  over   a  time   interval    sufficiently  short   that    long-term  effects 
may   be  neglected.      Since  the    averaging  period   is   not    standardized,    it 
must   be   specified;    it  varies   from  a   fraction  of   a  second  to   perhaps   one 
minute.      Long-term  stability   refers   to   slow  changes    in  average   frequency, 
and    is  usually   specified  for   crystal  oscillators   by   stating   fractional 
frequency  deviation  per  day,    e.g.  ,    5   parts    in   lOH   per  day.      Since 
atomic   frequency   standards   have   little   or   no   systematic  drift,    specifi- 
cations  for  them  will  often   be   given   in  terms  of   standard  deviation  or 
RMS  deviation  from  the  mean  over  a  specified  time   interval   such   as  one 
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day. 

Quartz  Crystal  Oscillators 

There  are  many  factors  which  affect  the  frequency  of  quartz  crystal 
oscillators,  which  must  be  taken  into  account  or  compensated  for.   In 
particular,  frequency  stability  is  a  function  of:  (a)  temperature;  (b) 
time  (aging);  (c)  stress,  vibration,  and  acceleration;  (d)  drive  level, 
and  (e)  nuclear  effects.   These  effects  are  discussed  by  Gerber  and 
Sykes  [29]. 

Ambient  temperature  changes  account  for  the  principal  frequency 
change  effects  in  most  crystal  units.   Current  practice  in  precision 
oscillators  is  the  use  of  a  double  proportional  oven  to  control  the 
temperature  of  the  crystal  unit.   This  provides  a  method  of  removing 
heat  from  around  the  crystal  at  the  same  rate  it  is  generated.   A  sim- 
pler, but  less  desirable,  method  uses  a  thermostatically-controlled 
oven.   In  this  method,  temperature  varies  between  a  lower  and  an  upper 
limit,  whereas  it  is  maintained  constant  when  using  proportional  con- 
trol. 

The  effect  of  time  (aging)  in  a  quartz  crystal  oscillator  is  easily 
observed  by  the  variation  in  frequency  deviation  from  day  to  day.   Accord- 
ing to  "Gerber  and  Sykes  [29]  ,  aging  of  thickness  shear  crystals  is  caused 
by  four  factors:  (a)  temperature  gradient  effects  lasting  several  min- 
utes to  several  hours  after  a  thermal  disturbance,  (b)  stress  relief 
effects,  (c )  change  of  mass  effects,  and  (d)  structural  changes  due  to 
imperfections  in  the  crystal  lattice.   In  all  crystal  oscillators  there 
seems  to  exist  an  initial  stabilization  period,  normally  lasting  a  few 
weeks,  in  which  the  oscillator  has  not  yet  attained  its  ultimate  fre- 
quency stability. 
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The  effects  on  frequency  stability  of  stress,  vibration,  accelera- 
tion, etc.,  may  be  reduced  by  careful  design  of  mounting  structures. 
Austin  [24]  has  made  an  experimental  study  of  the  effects  of  vibration, 
shock,  and  swinging  motion  on  frequency  stability  of  several  commercial 
crystal  oscillator  units.   He  concluded  that  stability  of  these  units 
was  degraded  negligibly  under  conditions  similar  to  those  which  might  be 
encountered  in  a  shipboard  environment. 

Another  factor  affecting  crystal  stability  is  drive  level,  or  ampli- 
tude of  vibration.  A  typical  2.5-MHz  crystal  unit  might  exhibit  a  change 
in  frequency  deviation  of  from  1  part  in  1CP  to  about  8  parts  in  10'  with 
a  change  in  drive  level  from  3  to  1000  microwatts  [29]  .  For  this  reason, 
automatic  gain  control  is  used  in  precision  crystal  oscillators. 

Nuclear  radiation  also  has  an  effect  on  the  frequency  of  crystal 
oscillators;  the  extent  of  this  effect  is  discussed  in  the  literature 
[13]  ,  [15]  ,  and  [26]  . 

Figure  5-1,  from  Gerber  and  Sykes  [29]  ,  shows  a  schematic  diagram 
of  a  modern  precision  crystal  oscillator.   The  double  oven  is  shown,  and 
proportional  control  would  be  used.   To  maintain  a  constant  drive  level, 
a  large  amount  of  negative  feedback  would  be  used  in  the  amplifier,  in 
addition  to  the  AGC. 

The  Sulzer  Model  D2.5  is  a  crystal  oscillator  of  the  type  shown  in 
Figure  5-1  and  was  in  use  at  the  Naval  Postgraduate  School  for  the  ex- 
perimental work  of  this  paper.   Its  drift  rate  was  determined  to  be 
about  2.5  parts  in  1011  per  day.   Other  units  of  this  general  type  often 
have  drift  rates  of  about  1  part  in  101*-  per  day,  while  similar  units 
used  to  stabilize  transmissions  of  some  U.S.  Navy  VLF  transmitters  have 
average  daily  drift  rates  as  low  as  a  few  parts  in  10    [29]  .   This 
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Figure  5-1.   Precision  Crystal  Oscillator 
general  type  of  oscillator  should  be  well  suited  for  shipboard  use 
in  a  VLF  relative  navigation  system. 

Atomic  Frequency  Standards 
It  is  interesting  at  this  point  to  note  the  difference  between 
a  primary  and  a  secondary  frequency  standard.   The  former  is  a  self- 
calibrating  device  which  is  intrinsically  accurate.   The  latter, 
however,  requires  calibration  against  a  known  precise  frequency 
before  it  can  be  considered  a  standard.   A  quartz  crystal  oscilla- 
tor may  thus  be  classified  only  as  a  secondary  frequency  standard. 
There  are  several  types  of  atomic  frequency  standards,  including  the 
hydrogen  maser,  the  cesium  beam  standard,  and  the  rubidium  vapor 
standard.   The  latter,  considered  a  secondary  frequency  standard 
since  its  frequency  initially  must  be  compared  with  some  other 
standard,  was  used  at  the  Naval  Postgraduate  School  as  an  adjunct 
to  the  preparation  of  this  paper.   The  hydrogen  maser  and  the  cesium 
beam  standard  are  both  primary  frequency  standards,   All  three  of 
these  atomic  standards  are  now  commercially  available. 
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Atomic  standards  in  general  derive  their  precise  frequencies  from 
frequencies  corresponding  to  the  energy  involved  in  the  transition  of 
electrons  between  energy  states  in  particular  types  of  atoms,  such  as 
hydrogen,  cesium  133,  and  rubidium  87.   The  cesium  and  rubidium  standards 
operate  on  the  same  basic  principles;  both  have  found  wide  use  in  VLF 
applications  as  well  as  other  uses.   Space  does  not  permit  detailed  des- 
criptions, but  since  the  rubidium  standard  was  used  in  some  of  the  ex- 
perimental work  for  this  paper,  its  basic  operation  will  be  described. 

The  heart  of  the  rubidium  standard  is  a  high-quality  quartz  oscilla- 
tor, whose  frequency  is  stabilized  by  a  passive  resonance  cell  filled 
with  rubidium  and  an  inert  buffer  gas.  The  operation  of  the  resonance 
cell  is  based  on  an  optical  pumping  principle.   The  rubidium  vapor  in 
the  cell  is  illuminated  by  filtered  light  from  a  rubidium  lamp.   A 
microwave  signal  is  applied  to  a  microwave  cavity  surrounding  the  reson- 
ance cell.   This  signal  is  obtained  by  multiplying  the  frequency  of  the 
quartz  oscillator.   When  the  frequency  corresponds  to  atomic  resonance, 
energy  transitions  are  induced  in  the  rubidium  atoms  and  a  photo  detec- 
tor observes  these  transitions  through  changes  in  the  amount  of  light 
absorbed  as  a  function  of  applied  microwave  frequency.   A  servo  loop 
connects  the  detector  output  and  oscillator,  locking  the  oscillator  to 
a  frequency  corresponding  to  the  center  of  the  resonance  line.   Since 
the  resonance  frequency  is  influenced  to  some  extent  by  the  pressure  of 
the  buffer  gas  in  the  resonance  cell,  this  pressure  must  be  properly 
adjusted  when  the  standard  is  manufactured.  The  cell  in  then  sealed 
and  the  unit  is  ready  to  serve  as  a  frequency  standard. 

Rubidium  vapor  frequency  standards  have  been  designed  for  U.S.  Army 
tactical  applications  and  for  guided  missiles  and  aircraft.   They  should 
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also  be   readily  adaptable  to   a   shipboard  environment.      State-of-the-art 
characteristics   for  these   standards,    according  to   McCoubrey    [J30]  ,    in- 
clude an  RMS  deviation  from  the  mean  of   5   parts    in   1012   for  a  one-day 
averaging  period   and  a  systematic  drift  of    less  than   3  X   10  per  month. 

Cesium  beam  frequency  standards    are  presently  being  used  as   the 
United   States   Frequency   Standard    (maintained  by  the  National   Bureau  of 
Standards,    Boulder,    Colorado)    and   to   stabilize  the  VLF   transmitters   of 
the  Omega  navigation   system  stations.      While   the   short-term  stability 
of  cesium  beam  standards    is   not   as   good  as  that  of  rubidium  standards, 
their   principal    advantages   are   excellent    long-term  stability  and  high 

intrinsic   reproducibility.      McCoubrey    [30]    gives   the   state-of-the-art 

1  ^ 
long-term  stability   as   2   X  10"    °  per  day    (RMS  deviation  from  the  mean). 

Stations  of  the  Omega  navigation  system  effect   an   intercomparison  of 

their  cesium  beam  standards    periodically    (currently   every   four   hours), 

and   are  thus   able  to   control   their   stability  to   better  than   1  X   10" 

per  day    [32]  . 

The  Varian  Model  V-4700    is   the   rubidium  vapor   standard   that  was 

used  at  the  Naval   Postgraduate   School   for   the  authors'    experimental 

work.      Over   a  period  of   about   four  months,    this   unit   exhibited   an  RMS 

deviation  from  the  mean  of   1.36  X   10"1    ,    using  an  averaging   period  of 

one  day. 
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6.   Antennas 

As  mentioned  in  Section  3,  the  long  path  signal  at  times  can  intro- 
duce serious  error  in  the  VLF  relative  navigation  system.   For  this 
reason,  it  is  desirable  to  use  some  means  of  discriminating  against  this 
component  of  a  received  VLF  signal.   The  cardioid  antenna  system,  which 
has  been  widely  used  in  radio  direction  finding  and  radio  compass  appli- 
cations, provides  a  ready  solution  to  the  problem.   The  basic  theory  of 
this  system  has  been  described  by  Terman  t3l  and  others,  and  will  be  re- 
iterated and  elaborated  upon  in  this  section.   The  loop  antenna,  basic 
to  a  cardioid  set-up,  will  first  be  discussed,  then  the  cardioid  system 
itself,  and  finally  the  results  of  measurements  on  antennas  at  the  Naval 
Postgraduate  School  will  be  presented. 

The  Loop  Antenna 

The  loop  antenna,  in  various  forms,  has  been  used  for  many  years. 
One  of  its  earliest  applications  was  reception  of  stations  in  the  broad- 
cast band  when  this  service  was  in  its  infancy.   Loops  used  for  VLF 
work,  as  well  as  most  other  applications,  have  a  diameter  which  is  small 
compared  to  a  wavelength.   A  loop  may  be  circular,  rectangular,  tri- 
angular, diamond -shaped,  or  some  other  similar  configuration.   Kraus  [2] 
states  that  far  fields  of  loop  antennas  of  different  shapes  are  the  same 
if  the  areas  enclosed  by  the  loops  are  the  same. 

A  vertical  loop  antenna  produces  a  polar  radiation  pattern  in  the 
familiar  figure-eight  shape,  as  shown  below  in  Figure  6-1.  The  xy-plane 
is  parallel  to  the  surface  of  the  earth.   Maximum  response,  of  course, 
occurs  when  the  plane  of  the  loop  is  oriented  parallel  to  the  direction 
of  arrival  of  the  electromagnetic  wave. 
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Figure   6-1.      Loop  Radiation  Pattern 

Appendix  III  describes   the  theory  of  operation  of   the   loop   antenna. 
Equations   are  derived  there   for   the   relative  response   and   phase   shift   as 
the  antenna   is   rotated   from  0°   clockwise  to    360°.      Figures    6-2   and   6-3 
show  plots   of   these  equations   for   0°   to   180°.      The  frequency   is   taken 
to   be   20  kHz,    and  the    loop   is   assumed   square,   with  the    length  of   a   side 
equal  to    1.635  m.       (This    length  was   chosen  because   it  was   the  measured 
length  of   a   loop  available   for   experimental   work.)     The    loop   is   assumed 
parallel    to  the  y-axis   of  Figure   6-1   and  direction  of   arrival    is   from 
right   to    left   along  the  x-axis    as   shown.      The  angle  of   rotation    0C    is 
taken   as   the   angle   between  the  normal   to   the   loop  and  the  positive  x-axis 
(this    is    shown   in  Figure   III-2) ,    so   that    0C    =  0°   and   180°   corresponds   to 
zero   response,   while    01  =   90°   and   270°   are  the  points   of  maximum  response. 

The  graph  of  Figure   6-2   was   arbitrarily  terminated   at   -35  db   for 
convenience,   since  output   theoretically  goes   to   -  oo  db   at    Cc  =   0°   and 
180°.      For   rotation   from  180°   to    360°,    the   amplitude   pattern   is   exactly 
the   same   as   for  the   first    180°   of   rotation   (this    is  obvious  from   in- 
spection of   Figure   6-1). 

For  the  phase   pattern,    the   second   180°  was   not   plotted   in  order   to 
show  greater  detail.      At    0C  =    180°   and    0C    =    360°,   there    is   a  discon- 
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Figure   6-3.      Loop   Relative   Phase   Shift 


49 


tinuity  of    180°.      The   pattern  for    0C    =   180°   to    (X    =    360°   looks   the  same 
as  Figure    6-3,    except    it   appears    inverted,    i.e.,    at    OL    =   180°,       C(?  — »» 
-180°;    at    OL   =   -270°,    t#    =   -(180°   +    3.4  x  10"4  deg);    and   at    a    =    360°, 
(/?—*-    -180°.      (The  symbol    — *-   is  used  to    imply  the  existence  of  the  dis- 
continuities  at    a    =   180°   and    tt    =    360°.)      The   180°   phase  reversal   on 
the   back  side  of  the   loop  antenna   is   significant,    for  one  must   be  care- 
ful  not  to   rotate  the   loop  antenna  to  the  wrong   side  of  the  null   if 
phase  measurements   are    important. 

The  Cardioid   Antenna  System 
A  cardioid  antenna  system  consists  of  a  loop  antenna  and  a  whip 
(or   "sense")    antenna  whose  outputs   are   appropriately  combined  to   pro- 
duce a  cardio id-shaped   radiation  pattern.      Figure   6-4   is   a  simplified 
schematic  of   such  a  system.     Both  the   loop  diameter  and   the  height  of 
the  whip  are   small  with  respect  to  a  wavelength.      Some  provisions   for 
tuning,   which  are  not   shown,   would  normally  be   installed  at   the  base  of 
the   loop  and   at   the  base  of  the  whip.      There  would   also   be   provisions 
for  varying  the  amplitude  of  the  loop  signal  or  the  whip  signal,   or  both, 
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Figure   6-4.      Cardioid  System  Schematic 
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as  well  as  provisions  for  varying  the  whip  phase  a  small  amount  above  and 
below  the  90°  phase  shift  shown.   These  provisions  are  necessary  to  allow 
adjustment  to  obtain  the  cardioid  pattern.   The  adjustment  consists  of 
making  the  whip  antenna  component  of  the  VLF  signal  equal  in  magnitude 
and  phase  to  the  loop  contribution,  for  the  case  in  which  the  plane  of 
the  loop  is  oriented  toward  the  distant  transmitting  station. 

Appendix  IV  contains  the  equations  necessary  to  plot  the  phase  and 
magnitude  of  the  cardioid  antenna  system  as  the  loop  is  rotated  about  a 
vertical  axis.   The  necessary  phasor  diagrams  are  shown,  and  Figure  IV- 3 
shows  the  phasor  diagram  for  proper  cardioid  adjustment,  with  the  loop 
oriented  toward  the  station.   Proper  operation  of  the  cardioid  system  is 
independent  of  the  location  of  the  whip  with  respect  to  the  loop,  so  long 
as  it  is  not  an  appreciable  fraction  of  a  wavelength  away. 

The  radiation  pattern  of  a  properly-adjusted  cardioid  antenna  sys- 
tem is  sketched  in  Figure  6-5.   Maximum  response  occurs  with  the  plane 
of  the  loop  antenna  oriented  parallel  to  the  x-axis. 


DIRECTION  OF 
SIGNAL  A&RIVAL 


Figure  6-5.   Cardioid  Radiation  Pattern 
Figures  6-6  and  6-7  are  plots  of  the  radiation  pattern  in  Car- 
tesian coordinates,  using  equations  derived  in  Appendix  IV.   The  plot 
is  for  a  frequency  of  20  kHz  and  a  square  loop  1.635  m  on  a  side.   Ro- 
tation angle  (X  is  the  angle  between  the  normal  to  the  loop  and  the 
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Figure  6-7.   Cardioid  Relative  Output  (-180  <0l<0  ) 
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positive  x-axis.      The  maximum  response   is   6  db,   which  merely   indicates 
that   the  voltage  delivered  to  the  receiver   is   twice  that   from  the   loop 
or  whip  alone.      The  minimum  response,   with  the   loop  pointed   away  from 
the  direction  of   signal    arrival,    is    approximately   -63  db.      In  the   forward 
direction,    6  db  beamwidth   is    180°,   providing  the  desirable  characteristic 
of   a  broad   front    lobe. 

It   should  be  noted   at   this   point   that   -63  db   is   not   the  minimum 
response  obtainable  on  the  "back"   side  of   a  cardioid   antenna   system. 
The   ideal  minimum  response   approaches   -oOdb  when  the  system  is  adjusted 
for  this.      The   adjustment  required   for  this   entails   adjusting  the  whip 
phase   for  minimum   system  response,    which  corresponds   to   setting  the  angle 
of  Es  to   +  <pa     rather  than  -  <p0     (see   Figures   IV-4  and    IV-7).      The   reason 
for   the  present   approach  to  the  adjustment   problem  is  that    for  a  weak 
station  it    is   nearly   impossible   accurately  to   adjust   the  whip  for  a  good 
null  while  the  antenna   is   pointed    in  the   reverse  direction;    this  has   been 
verified   in  the  laboratory. 

Phase   shift  of  the  cardioid   antenna  system   is   plotted    in  Figures 
6-8  and   6-9,    again  using  equations   derived   in  Appendix  IV.      It   should 
be  noted  that   the   T-axis   scales  of  the  two   figures   are  different.      For 
0°<&<  180°,   the  phase   shifts   no  more  than  0.3  x   10"2  deg   as   the   loop 

rotates.      At   a  frequency  of  20  kHz,    360°   phase  shift  corresponds   to   50 

9  -3 

ixsec,    so  0.3  x   10"^   deg  would  correspond   to    less   than   10       usee;    since 

normal   phase  variations  of  propagation  produce  changes  of    1   or  2  Usee 

on  a  phase   plot  of   a   received  VLF  station,   these   small  changes   would   not 

be   noticeable. 

An   interesting   aspect  of  Figure   6-9   is  the   phase  change  of   -90°    in 

the  vicinity  of    OC    =   -90    .      At    20  kHz,   this   would   correspond  to   a  phase 
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change  of  12.5  usee.      It  should  also  be  noted  that,  while  there  may  appear 
to  be  a  discontinuity  at  OC  =  0°  and  a  =   ±180°,  due  to  the  difference 
in  T-axis  scales  on  Figures  6-8  and  6-9  this  is  not  in  fact  the  case. 

Since  it  is  proposed  that  a  cardioid  antenna  system  be  used  aboard 
ship  in  the  VLF  relative  navigation  system,  there  has  been  some  concern 
that  the  sharp  phase  change  in  the  vicinity  of  the  cardioid  null  might 
cause  difficulties.   It  is  believed  by  the  authors  that  this  will  cause 
no  problems  for  the  following  reasons:  (a)  the  relative  broadness  of  the 
null,  and  (b)  experimental  evidence  that  the  90°  phase  shift  in  a  typical 
cardioid  system  does  not  occur  so  abruptly  as  the  theoretical  curve  of 
Figure  6-9  indicates.   The  latter  will  be  discussed  later  in  this  sec- 
tion.  The  former  is  seen  by  examining  Figure  6-7:  if  the  loop  is  rota- 
ted 15°  either  side  of  the  null,  cardioid  response  is  still  down  nearly 
35  db  from  that  of  the  forward  direction.   For  these  reasons,  it  is  felt 
that  no  special  measures  would  need  to  be  taken  aboard  ship  to  keep  the 
antenna  aligned  with  the  direction  of  the  transmitting  station,  other 
than  a  simple  servo  system  driven  by  the  ship's  gyro. 

Experimental  Work 

The  Tracor  Model  611  Cardioid  Unit  was  used  at  the  Naval  Postgraduate 
School.   This  unit  has  two  inputs:  one  directly  from  a  loop  antenna,  and 
the  other  from  a  whip  antenna.   Both  antennas  must  be  tuned  as  desired  in 
the  VLF  band  and  matched  to  50-ohm  coaxial  cable.   The  cardioid  unit  pro- 
vides adjustable  phase  shift  and  amplification  so  that  the  cardioid  pat- 
tern may  be  obtained  at  any  frequency  between  about  9  kHz  and  30  kHz. 
The  output  of  the  unit  is  switchable  between  cardioid,  loop  only,  and 
whip  only.   Although  the  cardioid  unit  is  specified  to  be  usable  with 
only  particular  types  of  Tracor  antennas,  it  was  found  to  be  usable  with 
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an   18-foot   whip  antenna  that   happened  to   be  on  hand   and  a  tuning  network 
designed   for  20  kHz.      The  complete  cardioid  system  manufactured  by  Tra- 
cor  comes   with   a  remotely -tuned   rotatable   loop  antenna  and   a  whip 
antenna  with  fixed  tuning.      This   system  was   used  to  obtain  data   pre- 
sented below. 

Radiation   patterns   for  the  cardioid   system  were  obtained   using  WVL 
at   20  kHz.      The   length  of   a  side  of   the  square   loop  was   1.635  m,    so   data 
is  directly  comparable  with  the  theoretical   patterns  of  Figures    6-6,    6-7, 
6-8,    and  6-9.      The   system  was   adjusted  by  setting  whip  amplitude  and  phase 
equal  to   that   of  the  loop  with  the   loop  oriented   toward   Ft.   Collins, 
Colorado   (the   location  of   WWVL) .      It  would  have  been   preferable  to   set 
up  a   nearby   low-powered   (milli-watt   or  micro-watt)    phase-stabilized 
transmitter  to   assure   steady   signal    levels  while  the   antenna  was   rotated, 
but   this   was    not   practicable. 

The  amplitude  data  obtained   is    shown  in  Figure   6-10.      Azimuth  of  0° 
corresponds    approximately  to  orientation  toward  the  transmitting  sta- 
tion.     The  numbers    in  parentheses   below  the  abscissa  correspond  to   angle 
of   rotation    (X  •      Due  to   difficulty   in  reading  the  rotator  unit,    relative 
azimuth  accuracy   is  probably   *3°.      Relative  amplitude   is   accurate  to 
better  than  ll  db.      The  two  major  sources  of   error  here  are  variations 
in  the  received  carrier   level  of  WWVL  and  calibration  of  the   receiver 
carrier   level   recorder  output. 

It   is    interesting  to   compare   Figure   6-10  with  Figures    6-6   and   6-7. 
Considering  that   real   antennas   never   perform  quite   so  well   as    ideal  ones, 
the   similarity   is   striking.      Overall,    the  null    is   somewhat   broader  than 
the  theoretical;    for  example,    the  output   is   35   db  down  at   about   22°   either 
side  of   the   null.      The   6  db   beamwidth   is    181°,   which   is  the  same   as   the 
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theoretical  value,  within  the  experimental  error.   The  two  points  appear- 
ing within  the  null  at  amplitudes  of  -32.5  db  and  -35.4  db  indicate  that 
the  cardioid  system  was  misad justed ,  with  whip  amplitude  slightly  great- 
er than  loop  amplitude.   A  careful  examination  of  Figure  IV-7  shows  that 
this  would  also  cause  phase  excursion  to  be  slightly  greater  than  >\/4. 
This  was  verified  experimentally,  as  seen  in  Figure  6-11. 

Figure  6-11  shows  phase  shift  data  obtained  as  the  loop  was  rotated, 
The  T  =  0  yUsec  line  is  an  arbitrary  reference  level.   Phase  data  is 
probably  accurate  to  il  or  2  Usee,  considering  accuracy  of  the  receiver 
and  propagation  changes  that  could  have  occurred  while  data  was  being 
taken.   The  smaller  figure  is  probably  acceptable,  since  data  was  taken 
during  the  daytime  when  the  ionosphere  was  stable,  over  a  period  of 
about  three  hours.   Angle  of  rotation  ot  is  again  indicated  in  parenthe- 
ses below  the  abscissa.   The  most  notable  difference  between  this  curve 
and  the  theoretical  data  of  Figures  6-8  and  6-9  is  the  broader  width  of 
the  90°  (  /\/4)  phase  change  in  the  vicinity  of  the  null,  as  noted  ear- 
lier.  The  rubidium  vapor  frequency  standard  was  used  as  the  local  phase 
reference  for  this  data. 

There  are  various  problems  associated  with  using  a  cardioid  antenna 
system  in  a  shipboard  environment.   One  is  the  problem  of  re-radiators 
located  near  the  loop-whip  system.   This  presents  a  problem  in  a  receiv- 
ing system  because  of  distortion  of  the  antenna  pattern,  principally  by 
resonant  lengths  of  metal  near  the  antenna.   Unfortunately,  it  was  not 
possible  to  investigate  this  problem  at  the  Naval  Postgraduate  School. 
It  is,  however,  the  opinion  of  the  authors  that  this  problem  would  not 
seriously  degrade  the  cardioid  antenna's  performance  in  the  VLF  relative 
navigation  system. 
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Another  problem  is  that  of  properly  adjusting  the  cardioid  unit. 
Whip  phase  and  amplitude  must  be  carefully  matched  to  that  of  the  loop 
at  each  frequency  used,  in  order  to  obtain  a  true  cardioid  pattern.   This 
may  be  difficult  for  distant  stations  whose  long  and  short  path  signals 
are  comparable.   Careful  attention  should  be  given  these  adjustments, 
since  depth  of  the  null  depends  directly  on  how  well  the  system  is  ad- 
justed. 
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7.   Determination  of  Signal  Strengths 

In  this  section  values  will  be  calculated  for  the  relative  strengths 
of  signals  received  via  the  long  and  short  great  circle  paths  as  a  func- 
tion of  the  time  of  day  and  day  of  the  year.   This,  along  with  the  phase 
of  both  signals,  is  necessary  in  order  to  calculate  the  phase  of  the  re- 
sultant signal. 

The  relative  strengths  of  the  two  signals  can  be  found  at  any  point 
on  earth  if  the  total  attenuation  on  each  path  is  known.   The  relative 
strength  of  the  two  signals  is  then  just  the  difference  of  the  short  and 
long  path  attenuations.   Upon  reexamination  of  the  electric  field  strength 
equation  in  Section  2  it  will  be  seen  that  an  assumption  is  being  made 
here  that  all  terms  except  the  last  two  in  this  equation  are  identical 
for  both  propagation  paths.   It  seems  only  reasonable  to  assume  that  most 
VLF  stations  have  a  symmetrical  radiation  pattern  and  that  the  amount  of 
excitation  in  two  opposite  directions  will  be  equal. 

Attenuation  rates  at  VLF,  mentioned  in  Section  2,  vary  from  about 
one  to  five  db/Mm  depending  on  the  time  of  day  and  propagation  path  con- 
ditions.  References  5,  8,  9,  18,  and  28  contain  many  of  the  curves  of 
attenuation  rates  versus  frequency  which  have  been  derived  from  experi- 
mental data.   These  curves  will  be  used  extensively  here. 

The  procedure  used  to  obtain  expected  attenuation  rates  at  any  par- 
ticular frequency  consists  simply  of  sifting  through  the  data  present  in 
these  references  and  extracting  the  attenuation  rates  given  for  different 
conditions  at  this  frequency.   From  this  information  estimates  of  the 
attenuation  rates  for  all  possible  path  conditions  and  times  of  day  can 
be  made. 

As  an  example,  at  13.6  kHz  the  following  values  are  given  in  the 
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references    stated. 
OC  (db/Mm) 


2.0 

1.8 
1.6 

3.5 

2.35 

1.95 

1.3 

1.0 

A*  =   1.0   to    1.5 

ACX  =  i.o 

&<*=    1.0 


conditions 


reference 

Taylor  [9]  pg.  1933 

Taylor  [9]  pg.  1933 

Chapman  and  Macario  C5] 
pg.  930 


sea,  night,  east  to  west   Taylor  and  Lange  [tQ 

pg.  609 


3.6 

land,  day,  east  to  west 

2.8 

sea,  day,  east  to  west 

2.6 

day 

day,    sea,   west  to   east 
night 


Taylor    [93    pg.    1933 

Chapman  and  Macario    £.5"] 
pg.    930 


sea,    day,   east   to  west  Watt   and  Croghan  [l8"2 

sea,   day,   north  to   south       Watt   and  Croghan   [18]] 


sea,    day,   west  to   east 
night,    sea 


Watt   and   Croghan  [l8D 
Watt   and   Croghan  C183 


sea,    night,    west   to  east,      Davies  C28D  Pg-    ^25 
15-20  kHz 


sea  to    land   change, 
15-20  kHz 

night   to   day  change, 
15-20  kHz 

W-E  to   E-W  change, 
15-20   kHz 


Davies  [28]  pg.  425 
Davies  f28]  pg.  425 
Davies    [28]   pg.    425 


From   these   attenuation  rates    (and   changes    in  attenuation  rates)    and 
other  discussions  of   attenuation  rates   such  as  those  by   Crombie  CO  an<^ 
Westfall  [2lj  reasonable   estimates  of   attenuation  rates    for    13.6  kHz   can 
be  made.      These  will   be   in  the  form  of  one   basic   attenuation  rate   plus 
additional   factors   which  depend  on  the  direction  of  propagation,    amount 
of  daylight   on  the  path  and   amount  of   land  on   the  path. 

The  values   chosen  are: 
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0<b   =1.0   db/Mm  basic   rate   (night,    sea,    W-E) 

A^n_(j   =    1.0   db/Mm  night   to  day  change 

^0<s_2   =  0-8  db/Mm  sea  to    land   change    (average   land) 

^^dN-S  =  °"1*  db/]ytn  W_E  to  N"s   direction   change 

^^dE-W  =   1-55  dVMm  W-E  to   E-W  direction  change 
The  change   in  attenuation  rate   as   a  function  of  direction  can  be 
approximated  by: 

A<XH(0)   =(  hMJt  )2 

d    ^  180° 

where    p    is    the  direction  of   propagation  with  respect  to  090°   true 
(0£  ^  £.  180°)  .      For  north  to   south   propagation   <J>    =   90°   and  for  east 
to  west   propagation    0    =   180°. 

The   attenuation  rate  over   any  path  can  be  calculated   from  the  fol- 
lowing  equation: 

*   =    *b   +   PD^"n-d    +  PL*°S-1    +^*V0> 
where  PD  =   fraction  of   the  path  which   is    in  daylight 

P-^  =   fraction  of   the  path  which   is  over   land. 

A  typical   plot  of  attenuation   for  a  24  hour  period    is    shown  in 
Figure  7-1.      The  path  here   is   Trinidad  to  Monterey    (short   path);    <fi   = 
145°,    P,    =  0.6   and  PD   is  obtained   from  the  computer   program  described    in 
Appendix   II.      The  total   attenuation   is   found  by  multiplying  0<   by  the 
path   length  obtained   from  the   geodesic   program    (Appendix  I).      The  cal- 
culated  attenuation  rates   vary  from  2.48  db/Mm  for  the  entire  path   in 
darkness   to    3.48  db/Mm  for  the  all-daylight   path. 

Figure  7-2   shows   a  comparison  of   the  long  and   short   path  attenua- 
tion from  Trinidad  to  Monterey.      The   long  path  here    is  considered  to   be 
in  the  opposite  direction  at   the  receiver   and  transmitter   from  the  short 
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Figure  7-1.  Calculated  attenuation  of  13. 6  kHz  signal  on 

Trinidad  to  Monterey  short  path  for  1  Oct.  1966 

Path  length  ■  665Q.Q5   km 
Kin.  =  16.72  db 
Max. -23. 37  db 
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Figure  7-2.  Calculated  attenuation  of  13.6  kHz  signal  on 
Trinidad  to  Monterey  long  and  short  paths  for 
1  Oct.  1966 

Short  path  length  =  6658.85  km 
Long  path  length  ■  33391.48  km 
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path  direction  of  propagation. 

It  should  be  noted  that  the  attenuation  rates  used  here  for  land 
paths  are  for  average  land.   Large  areas  of  very  high  or  very  low  conduc- 
tivity land  would  require  alteration  of  these  values.   Similar  care  must 
be  exercised  for  propagation  over  polar  ice  caps,  where  the  conductivity 
is  lower  than  that  of  sea  water. 

One  factor  not  mentioned  previously  is  the  loss  due  to  discontin- 
uities of  the  earth  conductivity  along  the  path  (Lj    in  the  electric 
field  strength  equation).  Westfall  L21J  mentions  this  factor  briefly 
as  having  a  magnitude  of  three  to  six  db  at  times.   Because  of  the  lack 
of  information  on  this  factor  it  has  not  been  taken  into  account  here. 

The  method  shown  in  this  section  could  be  used  at  any  other  fre- 
quency of  interest  in  the  VLF  range  to  obtain  similar  attenuation  rates. 
Some  complicating  factors  may  be  encountered,  however,  which  could  make 
the  calculations  more  difficult.   For  example,  at  16  kHz  it  was  found  that 
the  night  to  day  difference  over  land  is  greater  than  the  night  to  day 
difference  over  sea. 

Comparison  of  Experimental  and  Calculated  Values 

Figures  7-3  and  7-4  show  two  more  curves  of  calculated  attenuation 
for  the  Forestport  and  Haiku  Omega  stations  and  Figures  7-5  through  7-7 
show  the  observed  relative  carrier  levels  of  these  three  stations  as  ob- 
tained from  a  Tracor  Model  599G  VLF  phase  tracking  receiver.   No  measure- 
ments were  obtained  on  any  long  path  signals  or  on  the  Aldra  short  path 
signal  because  the  carrier  levels  were  below  the  level  of  control  of  the 
receiver  AGC,  giving  a  reading  of  zero  on  the  relative  carrier  level  out- 
put. 

Since  no  attempt  was  made  to  calculate  the  signal  strength  of  the 
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Figure  7-3.  Calculated  attenuation  of  13.6  kHz  signal  on 

Forestport  to  Monterey  short  path  for  1  Oct.  1966 
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Figure  7-5  •   Observed  relative  carrier  level  of  Trinidad 

13.6  kHz  at  Monterey,  30  Aug.  to  2  Sept.  19 66 

10  db/v>-  in.  vertical  scale 

night  average  -  day  average  =  4.04  db 


2    4    6    8    10   12   14   16   18   20   22   24 

UT 
Figure  7-6.   Observed  relative  carrier  level  of  Forestport 
13.6  kHz  signal  at  Monterey,  5  Sept.  to  8  Sept- 
1966. 

10  db/ln.  vertical  scale 
night  average  -  day  average  =  2.47  db 
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Figure  7-7.  Observed  relative  carrier  level  of  Haiku  13.6 

kHz  signal  at  Monterey,  2  Sept.  to  5  Sept.  1966 
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various  stations  and  this  information  is  not  readily  obtainable  from  the 

Tracor  receivers,  the  only  quantities  that  can  be  compared  meaningfully 

are  the  difference  between  the  carrier  levels  for  the  entire  path  in 

darkness  and  in  daylight  and  the  day-night  attenuation  difference.   While 

the  attenuation  decreases  at  night  and  the  carrier  level  increases  during 

the  same  period  the  difference  between  the  day  and  night  values  should 

be  equal  since  they  are  complementary  quantities.   (This  is  easily  seen 

in  the  electric  field  strength  equation  of  Section  2  where  &<d/106  is 

the  attenuation  and  E„  fc  is  the  carrier  level.) 

z  ,n 

A  comparison  of  the  calculated  and  observed  values  of  the  night -day 
difference  is  shown  below. 

Station  Calculated  difference    Observed  difference 

Trinidad  6.65  db 

Forestport  4.01  db 

Haiku  3.86  db 

The  values  listed  under  (a)  were  observed  with  the  receiver  servo 
time  constant  set  at  150  seconds  and  those  under  (b)  with  a  15  second 
time  constant.   The  lower  differences  observed  with  the  150  second  time 
constant,  corresponding  to  Figures  7-5  and  7-7  are  apparently  because 
of  the  erratic  nature  of  nighttime  phase. 

Calibration  of  the  signal  level  scale  used  on  the  recorder  was  very 
difficult,  leading  to  some  possible  inaccuracies  in  these  measurements. 
(Calibration  was  complicated  by  the  short  duty  cycle  of  the  transmitting 
stations;  10%.)   Thus  while  these  measurements  do  provide  some  measure  of 
confidence  in  the  calculated  attenuation  rates,  it  is  believed  that  more 
comprehensive  measurements  of  the  carrier  levels  of  these  signals  would 
be  necessary  if  an  accurate  check  of  these  attenuation  rates  was  desired. 
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8.   Calculation  of  VLF  Phase  Velocities 

Figures  2-8  through  2-11  show  several  curves  of  calculated  phase 
velocities  for  VLF  signals  for  different  conditions  which  could  be  used 
directly  to  obtain  the  desired  phase  velocities.   However,  it  would  be 
more  convenient  and  accurate  to  make  use  of  an  equation  for  phase  velocity 
which  takes  into  account  the  frequency  and  conditions  of  the  ionosphere 
and  the  earth.   Such  an  equation  is  available  in  the  Omega  System  book 
[3l],  page  3-20: 

V              h-                          v      1 2 
_JL_  =1-0.36  — i-  +   (2-ITn  -  (A      -  ip  ,)  ^— 

Here  V  =  phase  velocity  of  the  VLF  signal 

VQ  =  free  space  velocity  of  the  wave 

=  3  x  10  m/sec 

0.36   =   a  constant   chosen  to   give   a  good  fit    in  the   10-20   kHz 
frequency  range 

h-    =   ionospheric  height   for  which  the  conductivity    (tf")   =   10 
mho  s/m 

a  =   radius  of   the  earth  =   6378.155  km 

*p      =   phase  shift   upon  reflection  at  the  ground 

(D-    =   phase  shift   upon  reflection  at  the   ionosphere 

n  =  mode  number    (=    1   for  primary  mode) 

f   =   frequency 

Some  values  are  given  for  the  phase  shifts  that  can  be  expected; 

*~P^   =  -160°  for  daytime  propagation,  -165°  at  night,  -165°  for  daytime 

propagation  to  the  east,  and  -155°  for  daytime  propagation  to  the  west. 

Assuming  that  the  -160°  phase  shift  is  for  north-south  propagation,  which 

seems  reasonable,  the  values  for  nighttime  propagation  to  the  east  and 

west  would  then  be  -170°  and  -160°  respectively.  The  ground  phase  shift 
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is  zero  for  sea  water  and  the  values  given  for  land  are  -0.09  radians 
during  the  day  and  -0.13  radians  at  night,  for  a  ground  conductivity  of 
10"2  mhos/m. 

Using  the  given  ionospheric  heights  of  68  km  during  the  day  and  86 
km  at  night  (derived  from  an  exponential  conductivity  profile  and  taking 
the  height  for  which  KT   =  10~°  mhos/m)  the  phase  velocity  curves  of  Fig- 
ure 8-1  were  calculated.   A  comparison  of  these  curves  with  similar 
curves  in  Section  2  (Figures  2-9,  2-10,  2-11)  shows  a  very  close  corres- 
pondence.  Several  experimental  phase  velocities  from  Figure  2-14  are 
also  plotted  here  for  further  comparison. 

The  phase  velocity  at  any  point  along  a  path  is: 

V.  (x)  =    2*f 


P        dO/dx 

where    &    is   the  phase  in  radians  at   the  point   x  on  the  path.      The  total 
phase   along  the  path   is   then: 

.d 

d> 

:xj 


o  ■  2jrtf  -& 


o  'P 

If   the  phase  velocity   is   constant  this   becomes: 


e  = 


2TTfd 
VP 


Using  a  method  similar  to  that  employed  in  determining  the  path 
attenuations  of  Section  7,  an  average  phase  velocity  can  be  obtained 
for  each  path.   It  can  be  calculated  by  the  following  equation: 

p    pb      pd   ^     D4*  pN-D    S   pl-s 

where  V  b  =   basic   phase  velocity;    phase  velocity  for   night,    land, 

W-E  from  Figure  8-1. 

&V  j(0)   =    increase   in  phase  velocity  due  to  direction  of   propagation. 
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PD  =  proportion  of  path  in  daylight. 
AV  N   =  increase  in  phase  velocity  for  daylight  propagation. 

Pg  =  proportion  of  path  over  sea. 
AV  -,    =  increase  in  phase  velocity  for  propagation  over  sea. 

The  increase  in  phase  velocity  with  direction  can  be  obtained  from 
the  following  equation: 


AVpd(*)  =  AV^-jIb- 


where    0    is   the  direction  of   propagation,   measured   as    in  Section  7.      From 
Figure   8-1,    the  differences   between  the  phase  velocities   for   E-W  and  W-E 
propagation   are  different   for  daytime  and  nighttime  propagation.      £V  o 
for  daytime  can  be  approximated  by  the  average  daytime  differences    in 
phase  velocity  with  direction  and  the   nighttime  value  can  be  obtained 
similarly.      This  can  be  expressed  by  the  following  equations: 

AVpd   =  ^(VpsDEW  -   VpsDWE   +  VplDEW  "   VplDWE)   ^Y^™ 

AVpd   =  3sCVp8NEW  -   VpsNWE   +   VplNEW  -  VplNWE)nighttime 


Here     vdst)etj  =  phase  velocity  for   sea,    day,    E-W 

V  -ifirjp   =  phase  velocity  for   land,    night,   W-E  and   so  on. 

Using  reasoning   similar  to   the   above   it    is   found  that   two  values 
for   &V  N  n  must   be  used   (one  for   E-W  propagation  and  one   for  W-E)   and 
one  value  for    AVnl    0   can  be  used.      These  are  given  by  the   following 

p  J.  ~  s 

equations  : 

A  VpN-D   =  ^(VpsDWE   "   VpsNWE    +  VplDWE   "   VplNWE) '   W"E 
^VpN-D   =  ^(VpsDEW  "   VpsNEW  +   VplDEW  ~   VplNEW)  »    E_W 

^Vpl-s   =   ^(VpsDWE   "   VplDWE   +  VpsDEW  "  VplDEW 
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+   VpsNWE   "   VplNWE   +   VpsNEW        VplNEW) 

The  range  of  the  phase  velocity  changes  being  averaged  to  obtain 
AVn;L_s  i-s  about  4  x  10"  V _/V  .  The  range  of  phase  velocity  changes 
being  averaged  for  AV  N_D  and  AV  ^  is  smaller,  indicating  that  the 
error   in    AVD^_S    is   the    largest  that  would  be  encountered. 

The  curve   shown   in  Figure   8-2    is   a  typical    24  hour   plot   of  phase 
delay  calculated  by  the   above  method.      This   is   for  the  Trinidad  to  Mon- 
terey  short   path  with  the   same   parameters   as   used   in  Section   7(0  = 
145°,   P     =   0.4,   PD  from  computer  program  of  Appendix   II).      The   13.6  kHz 
phase  velocities  used  to   calculate  this   curve  are  listed  below. 

VpsDWE  =    I- 00002 13VQ 

VpsDEW=    1- 00<W27  3V0 

V   BKnn  =  0.9974367V 
psNWE  o 

VpsNEW=   0.9976842Vo 
VplDWE  =   0.9998199Vo 
VplDEW  =    l'0002154Vo 
VplNWE   =   0.9972606Vo 
VplNEW=0'9974986Vo 
Comparison  of   Experimental    and  Calculated  Phase   Delays 
Figures   8-3  through  8-6   show  both  calculated   and  observed   24  hour 
phase  delay   curves.      The   average  of  the  observed  daily  curves    is   shown 
shifted  by    1.5   inches    (15  or   30/<sec)    for  convenience  and   ease  of   com- 
parison with  the  calculated  curves.      The  total   phase  delay  of  the  signal 
from  transmitter  to   receiver   is   not   easily  obtained  experimentally;    thus 
only  the  day-to-day  changes    in  the   phase  delay  are  compared  here. 

Figures   8-5   and   8-6   show  a  very  close  correspondence   between  calcu- 
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lated  and  observed  curves.   It  is  interesting  to  note  that  there  appears 
to  be  a  time  lag  in  the  observed  phase  delay  during  periods  of  rapid 
change,  which  could  be  due  to  the  long  time  constant  used  in  the  receiver 
(150  sec).   The  lag  also  seems  more  pronounced  on  stations  with  lower 
signal  strengths  at  the  receiver. 

Figure  8-4  shows  differences  in  the  calculated  and  observed  curves 
which  seem  to  point  towards  some  error  in  the  calculations.  The  ampli- 
tude and  shape  of  the  observed  diurnal  shift  are  believed  to  be  approx- 
imately correct. 

Figure  8-3  is  a  good  example  of  a  diurnal  shift  that  is  not  a  trap- 
ezium.  Propagation  over  a  polar  path, as  is  the  case  here,  yields 
diurnal  shifts  similar  to  those  shown  in  Figure  8-7,  and  are  markedly 
different  during  different  times  of  the  year.   The  diurnal  shift  on  polar 
paths  is  discussed  more  thoroughly  by  Blackband  Cl2].   In  Figure  8-3  the 
poor  correlation  of  the  calculated  and  observed  curves  may  be  partially 
due  to  the  weak  signal  strength  of  the  Aldra  station  at  Monterey.   The 
phase-tracking  receiver  did  not  seem  to  "lock  on"  to  the  received  phase 
until  a  large  change  in  phase  delay  occurred,  such  as  encountered  during 
sunset  along  the  path.   This  can  be  seen  during  the  initial  portion  of 
the  readings,  started  at  1900  on  23  September. 

Figure  8-8  shows  a  comparison  of  an  experimental  phase  velocity 
curve,  obtained  by  Thompson  et  al  [lb"]    in  Sydney,  Australia  at  16  kHz, 
with  a  calculated  curve.  The  data  used  for  the  calculated  curve  is 
approximate  (position  of  the  receiver,  proportion  of  land  on  the  path, 
etc.)  and  in  view  of  this  the  close  correspondence  of  these  two  curves 
is  remarkable. 
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Figure  8-2*  Calculated  phase  delay  on  Trinidad  to  Monterey 
short  path  at  13.6  kHz- 
Minimum  value  is  22195 . 759/^sec  delay 
10y*sec/inch  vertical  scale 
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Figure  8-3-.  Aldra  to  Monterey  phase  delay,  short  path 
13.6  kHz,  23-26  September,  1966 
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Elgure  8*-4.  Haiku  to  Monterey  phase  delay,  short  path 
13.6  kHz,  25  July-1  August,  1966. 
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Figure  8-5,  Forestport  to  Monterey  phase  delay,  shoi^  path 
13.6  kHz,  6-9  October,  1966 
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Figure  8-6.  Trinidad  to  Monterey  phase  delay,  short  path 
13.6  kHz,  11-14  October,  1966 
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Figure  8-7.  Typical  diurnal  shifts  on  a  polar  path. 
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Figure  8-8.  Observed  and  calculated  phase  of  GBR  at  Sydney, 
Australia  for  October  1961  and  October  1966 
respectively. 
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9.   Calculation  of  Resultant  Signal  Phase 

As  discussed  in  Section  3,  the  amplitude  and  phase  of  a  signal 
consisting  of  long  and  short  path  components  can  be  obtained  by  the 
addition  of  two  phasors.  In  the  navigation  system  proposed  in  this 
paper  the  phase  of  the  resultant  signal  at  a  receiver  is  of  primary 
interest.  The  amplitude  of  the  signal  is  of  interest  only  insofar  as 
it  is  sufficient  to  allow  phase  tracking. 

The  phase  of  the  resultant  signal  can  be  found  from  the  phase  of 
the  long  and  short  path  signals  and  their  relative  amplitudes  or  atten- 
uations.   First  take  the  magnitude  of  the  phasor  E..  in  Figure  3-3  to  be 
unity,  or  zero  db.   Then  the  magnitude  of  E~  is  the  relative  attenuation 
of  the  short  path  signal  with  respect  to  the  long  path  signal.   It  can 
reasonably  be  assumed  that  when  )EJ<  0*1  the  phase  of  E,  is  not  appre- 
ciably affected  by  E.  (for  ff2   =  0.1 /9Q°.  <p   =  5.7°  in  Figure  3-3). 

The  relative  attenuation  for  this  case  is  given  by  the  following 

equation: 

|E2) 
A  =  20  log  -I-£j-  =  -  20  db. 

Thus  the  long  path  signal  must  be  taken  into  account  whenever  the 
strengths  of  the  two  signals  are  within  approximately  20  db  of  each 
other. 

To  obtain  the  phase  of  the  resultant  signal,  the  phase  delay  in 
microseconds,  calculated  in  Section  8,  must  be  converted  to  an  angle. 
This  is  easily  done  by  converting  to  wavelengths,  where  one  wavelength 
equals  1/f  seconds  (at  20  kHz,  1/f  =  50/<sec).   One  wavelength  corres- 
ponds to  360°  and  6  is  then  the  difference  between  the  phase  delays  on 
the  long  and  short  paths.  The  angle  0  is  then  given  by  the  following 
equation: 
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Eo   sin  9 
tan  0  = 


1  +  |E2|  cos  0 

This  angle  could  also  be  found  by  adding  graphically  the  phasors  E^  and 
E2. 

It  must  be  remembered  here  that  the  total  phase  delay  of  the  re- 
sultant signal  is  equal  to  0  plus  the  phase  delay  of  the  short  path  sig- 
nal . 

Figures  9-1  through  9-5  show  the  calculated  curves  of  long  and  short 
path  attenuation,  and  short  path  and  resultant  phase  at  Sydney,  Australia, 
and  Wellington,  New  Zealand.   The  transmitter  station  used  was  GBR  at  16 
kHz  and  the  day  was  assumed  to  be  2  3  September  1966. 

In  Figure  9-2  the  largest  difference  between  the  short  path  phase 
and  the  resultant  phase  is  slightly  less  than  one  microsecond  and  the 
difference  in  attenuations  at  that  time  is  20  db.   This  verifies  the 
earlier  assumption  that  the  resultant  phase  differs  from  the  short  path 
phase  (by  an  appreciable  amount)  only  when  the  strength  of  the  long  path 
signal  is  within  20  db  of  the  short  path  signal. 

Figure  9-4  shows  a  case  where  the  resultant  phase  is  considerably 
different  from  the  short  path  phase.  The  largest  difference  here  is  12.1 
/Asec  (at  0300)  or  about  70  electrical  degrees,  and  does  not  occur  at  the 
time  when  the  long  path  signal  is  strongest  (long  path  signal  is  0.6  db 
weaker  at  0  300,  1.8  db  stronger  at  0500).   From  about  0800  to  1800  there 
is  approximately  20  db  of  separation  between  the  two  signals  and  the  re- 
sultant phase  is  very  close  to  the  short  path  phase  as  expected. 

In  Figure  9-5  the  latitude  of  Wellington  was  changed  by  four  min- 
utes, lengthening  the  short  path  by  about  seven  km  and  shortening  the 
long  path  by  the  same  amount.   The  total  phase  delay  on  the  short  path 
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Figure  ^-Iv  Attenuation  on  GBR  to  Sydney  short  and  long  paths 
at  16  kHz. 


87 


.    .              1 

/        \ 

1 
1 

1 

1 
I 

1 

1 

V 

/ 

r 

\ 
\ 

rat 
/ 

\ 

) 

r. 

/ 
/ 

A 

T 
1 

20/Asec 

\ 

Phase 

i 
/ 

/    l        y 

i 

V 

■\ 

/ 

i 

i 

A A 

Pi        n 

i 

! 
i 
resultant 

i 

U       -Id 

sho 

rt  pa^n 

i 

&r~f-3 re      c 

± _ —    . — i 

i 
1 

8 


12 


UT 


16 


20 


24 


Figure  9-2.  Short  path  and  resultant  phase  of  GBR  at  Sydney 
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Figure  9-&.  Short  path  and  resultant  phase  of  GBR  at  Wellington 


90 


Phase 


Figure  9-5 .<  Short  path  and  resultant  phase  of  GBR  at 
Wellington.  (Short  path  7  km  longer) 
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is  increased  slightly  as  a  result  of  this,  but  the  diurnal  shift  remains 
the  same.   The  resultant  phase,  however,  is  considerably  different,  show- 
ing that  a  change  in  position  has  a  larger  effect  on  the  resultant  phase 
than  on  either  the  long  or  short  path  phase  alone. 

No  experimental  data  has  been  found  with  which  these  calculated 
curves  could  be  directly  compared;  however,  Thompson  et  al.^14}  and 
Crombie  C^D  have  published  articles  which  discuss  this  phenomenon  and 
include  some  experimental  evidence  taken  at  Sydney  and  Wellington. 

It  should  be  noted  that  the  resultant  phase  curves  as  shown  here 
could  not  be  observed  by  using  a  loop  antenna.   The  180°  phase  reversal 
with  a  reversal  of  the  direction  of  the  incoming  signal  would  cause 
some  different  resultant  phase  to  appear.   This  could  be  taken  into 
account  in  the  calculated  curves,  however,  by  adding  180°  to  the  phase 
of  the  long  path  signal  before  adding  the  phasors. 

The  concept  of  phasors  as  discussed  here  provides  another  means  of 
describing  the  phenomenon  of  cycle-slipping.   Figure  9-6  shows  the  am- 
plitude and  phase  of  some  station  over  a  24-hour  period  as  represented 
by  a  phasor  trace. 


Figure  9-6.   Phasor  representation  of  a  24-hour 
phase  and  amplitude  curve 
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This  figure  is  drawn  assuming  a  trapezoidal  diurnal  pattern  of  both 
phase  delay  and  amplitude,  both  being  greater  at  night.   Thus  A  repre- 
sents  the  daytime  phase  and  amplitude,  B  represents  the  nighttime  phase 
and  amplitude,  8"  is  the  diurnal  phase  shift,  and  |B|  -  |A|  is  the  night- 
time increase  in  amplitude. 

Cycle-slipping  has  occurred  if  the  curve  traced  by  the  phasor  en- 
closes the  origin,  as  shown  by  the  dotted  line  in  Figure  9-6.   This 
could  be  due  to  long  path  interference,  loss  of  signal  (in  which  case 
the  magnitude  of  the  vector  is  zero),  or  several  other  causes.  The 
case  for  which  the  diurnal  shift  is  greater  than  one  cycle  (2 TO  should 
not  be  confused  with  cycle-slipping  as  shown  here.   Figure  9-7  shows  one 
possible  curve  of  this  type  and  it  can  be  seen  here  how  easily  cycle- 
slipping  could  occur. 


Figure  9-7.   Phasor  representation  of  a  24  hour  phase  and 

amplitude  curve  with  a  large  diurnal  phase  shift 

A  rapid  change  in  phase,  as  sometimes  experienced  during  sunset  or 
sunrise,  could  cause  the  carrier  level  in  the  receiver  to  be  greatly  re- 
duced for  a  time,  particularly  if  a  long  servo  time  constant  is  used. 
This  could  result  in  a  phasor  trace  similar  to  the  one  shown  in  Figure 
9-8,  with  an  increased  probability  of  cycle-slipping.   This  rapid  phase 
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change  and   resultant  cycle-slipping   is  typical  of   north-south  propagation 
paths. 


Figure   9-8.      Phasor   representation  of   loss  of  carrier 
level   due  to  rapid  phase  change 
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10.      Conclusions 

Attenuation  rates   and  relative   signal   strengths  of    long  and   short 
path  signals   can  be   determined   easily.      The  accuracy  of  the  method   used 
here  has    not   been  determined,   however,   because  of  a   lack  of   dependable 
experimental   data. 

It  has  also  been  demonstrated  that  the  phase  of  a  signal  can  be 
predicted  for  a  typical  24  hour  period,  both  with  and  without  long  path 
interference.  Experimental  data  has  shown  the  short  path  phase  calcula- 
tions to  be  very  accurate  for  some  cases  when  three-to  five-day  averages 
of  observed  phase  are  considered.  Lack  of  experimental  data  prevented  a 
check  of  the  calculated  phase  in  the  presence  of  long  path  interference; 
however,  it  was  shown  that  when  this  type  of  signal  is  present  the  phase 
is  highly  sensitive  to   a  change   in  receiver  position. 

Additional   experimental   data,    possibly   from  some   location  other 
than  Monterey,   California,   would   be   highly  desirable  and  would  probably 
lead  to   some  modifications  of  the  parameters   used   in  the  calculations 
presented  here.      It   should  be  emphasized  that   these  calculations   are  a 
first    iteration   in  an  attempt' to  predict,    at   any  point  on  earth,    the 
phase  of  a  VLF  signal. 

A  card io id   antenna   system   should  be   suitable   for   use  aboard   ship 
in  a  VLF  relative  navigation  system.      A  simple   servo   system  driven  by 
the   ship's   gyro,    and   provided   at    intervals   with  the  true  bearing  to   the 
desired  VLF  station,    should  be   adequate  to   keep  the  cardioid  pattern 
oriented  with   sufficient   accuracy.      With   proper   adjustment  of   the  car- 
dioid  system,    interference   from  the   long  path  signal  may  be  eliminated. 
In  case  of   a   location  at   which  the   long  path  signal    is   always   equal   to 
or   greater   in   strength  than  that  of   the  short   path,    however,    it   may  be 
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impossible  to  adjust  the  cardioid  system  properly,  and  hence  it  could 
not  provide  the  desired  amount  of  long  path  signal  attenuation. 
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APPENDIX   I 
GEODESICS    AND   LONG  PATHS 
A.      Geodesies 

The  great   circle  distance  between  two   points   on  a  sphere   is   easily 
calculated;    however,    determination  of  the  length  of   a  geodesic    (analo- 
gous  to   a  great   circle)   on  an  ellipsoid   is  more  difficult. 

Kirkland   C2CQ   states   very   concisely  the  simplest  method  of  obtaining 
the   length  of   a  geodesic. 

Given  two   points    (&i>    A-,),    (#2'    *  ~)   find   their  respective  para- 
metric  latitudes,  j$  -^   an^   B  2   as   shown   in  Figure   1-1. 

0  =  geodetic   latitude 
/3  =  parametric   latitude 

0.  =  equatorial   radius 

v  =  polar  radius 


tan  0    =   -  . 

z  =   a   cos  A 

y  =  b   sin/? 

dz     _   _      a 
dy     ~     "     b 


dz 
dy 


tan^? 


tan  B      -  tan  0  =   c   tan^J 

sec^^   =   1   +  c2   tan2  0 
1 


Figure   1-1.      Geodetic   and 
parametric   latitudes 
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The  arc   distance  between  the  two   points    (x)    is   then   found   from: 


cos   x   =   sin 
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fi  2   cos    (<\l   -  A2) 
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and    the   great   circle  distance   is  D  =   a  x   . 

Now  using  the  correction  factors   derived  by  Lambert  CO  the   length 
of  the  geodesic    is: 

d    =   D   -   AQP  -   BQQ 

(a  -   b)    (x  -   sin  x) 


where  P  = 


Q  = 


4(1   +  cos   x) 
(a  -  b)    (x  +   sin  x) 


4(1   -  cos  x) 
AQ   =    (sin  /2  l   +   sin  p  2)2 

BQ   =    (sin^  !   -   sin  ($  2)2 

The   accuracy  of  this  method   is,    according  to   Lambert,    better  than 
*  0.05   km  but   breaks   down  near   the   antipodes    (i.e.    for  x  >175°).      This 
method  was   checked   against    the  examples    in  an  article   by  Rainsford  C4] 
with  the  largest  error  being  about  0.85  km.      (It   should  be  noted  that 
Rainsford   specified  the  position  of  the  points  to   within  0.0001   second.) 
B.      Long  Path  Distances 

The  method  used   in  this   paper  to  obtain  the   long  path  distance  be- 
tween two   points    is  to   find   the   length  of   the   perimeter  of  the  ellipse 
on  which  the  two   points    lie    (analogous   to   the  great   circle  on  a   sphere) 
and  subtract   the   geodesic   length  from  this  distance. 

The  major  axis  of  this   ellipse   is   the  equatorial   radius  of  the 
ellipsoid   and  the  minor  axis   can  be   found   as   follows: 

/3   =  parametric   latitude 
x     =  arc  distance  between  the  two   points 
cos  Dt    =  -£°s  A  ""  coS  x  COS   C 


sin  x  sin  C 
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Pole 


in  /Q 


-  cos   x  sin 


fi 


sin  x  cos  Q  -, 


Equator 

Figure  1-2.      Ellipse  of  geodesic 


cos    #     =   sin 


(X     COS   ^  ,     =    COS  A     i  -1/  1    -       <— 4 7! Cl^U 

K  1  ^    [   V  sin  x  cos  fi  Y         J 


X    COS   & 

Here  gamma   is    the  parametric   latitude  of  the  northernmost   point   on  the 
ellipse.      From  Figure  1-3: 

z  =  a  cos  y 

z2  y2        . 

77-    +  ~        ~     1 


2  I2 

y       =    b 


ll_   Z2    =    b2(l    -    COS2    tf    ) 


Figure   1-3.      Parametric   and 
geocentric    latitudes 


The  minor   axis  of   the   ellipse  of    interest   is   then: 


=  -\/z2   +  y2'    =    l/b2    +   << 


V   -  b*)    cosz)J 

The   length  of  the  perimeter  of   the  ellipse  described  by  a  and   r   is 
found  by  evaluating  the  elliptic    integral: 


E(k)    =  J      -\fi   -  k2   sin2x"   dx 


where 


=  VTTE^ 
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The  perimeter  is  then: 

P  =  4  a  E  (k) 
and  the  long  path  distance  is: 


dL  =  P  -  d 


where  d  is  the  geodesic  length. 

The  FORTRAN  60  program  used  to  compute  these  distances  is  shown  on 
the  following  pages  along  with  sample  results. 
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APPENDIX  II 
PROPORTION  OF  DAYLIGHT  AND  DARKNESS  ON  PATHS  OF  PROPAGATION 

The  difference  in  the  height  of  the  D-layer  of  the  ionosphere  from 
day  to  night  is  the  largest  single  cause  of  variations  in  the  phase  vel- 
ocities of  VLF  waves  and  causes  variations  in  the  attenuation  rates  of 
a  magnitude  equal  to  those  caused  by  other  factors.   Because  of  this  the 
percentage  of  daylight  and  darkness  on  each  path  of  interest  must  be 
calculated  as  a  function  of  time. 

This  problem  has  been  solved  by  Brady  and  Crombie  C22^.   However, 
since  the  following  solution  was  developed  independently,  it  will  be  out- 
lined here. 

The  following  is  based  on  the  premise  that  the  percentage  of  day- 
light on  a  path  can  be  found  if  four  points  are  known:  (1)  transmitter 
location,  (2)  receiver  location,  (3)  and  (4)  intersection  of  the  sunrise 
and  sunset  lines  with  the  great  circle  path  between  (1)  and  (2).   (A 
spherical  earth  of  unspecified  radius  is  assumed.) 

The  great  circle  path  can  be  described  by  the  longitude  of  one  of 
the  points  where  it  crosses  the  equator  (Ai)  and  its  declination  at 
the  equator  (0<t).   The  sunrise-sunset  line  forms  another  great  circle 
that  can  be  described  in  the  same  way.   The  intersection  of  these  two 
great  circles  on  the  sunrise  side  of  the  earth  will  determine  the  coor- 
dinates of  point  (3)  mentioned  previously.   The  geometry  of  the  situation 
is  shown  in  Figure  II-l. 

Here: 

\      =  15t  -  90  degrees 

OC  2  =  °*s  +  90  degrees 
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where 


t  =  GMT 
0<.   =  declination  of  the  sun 


sunrise  line 
:ircle  path 


\3        n*  "i 

Figure   II- 1.      Point   of   sunrise  on 
great   circle   path 

Using   spherical    trigonometry  the  coordinates  of  P   (^-,    \  0    can 

easily  be   found: 

cos^Cot!    -  DC2) 


a  =   tan 


1     [tan1- 


(A2  -  Xi> 


cos^C^!    -  DC2) 


+  tan 


1     tan^(A2   -\J      —       V- 

I  z         M      sin%(X  !   +  X2> 


tan   c   =  tan  a  cos(180°   -   &C„)    =   -  tan  a  cos  &<  2 


+   c 


A3   -    A  2 

n  ^  3   =   sin   a   sin  0^2 


si 

The  point  at  which  the  sunset  line  crosses  the  great  circle  path 
is  then: 

A  4  =  X  3  -  180° 


£4  --03 


At  a  height  of  70  km  the  point  P  will  be  further  along  the  great 
circle  path  in  a  direction  away  from  the  sun.   Figure  II-2  shows  the 
arc  distance  by  which  the  sunrise  and  sunset  lines  are  shifted. 

Here: 

h  =  ionospheric  height 


110 


S  =  screening  height 
R  =   earth  radius 


sun's   rays 
Figure   II-2.      Sunrise/sunset   at    ionospheric  height 
The   screening  height,    mentioned  by  Brady  and   Crombie  C223»    is  the 
height    below  which  the  components  of  the   sun's   rays  that  most   affect 
ionization  cannot   penetrate.      This,    in  effect,   changes   the  radius  of 
the  earth  by  an  amount   S.      From  Figure  II-2  then: 


cos 


f  = 


R  +  S 


R  +  H 

The  point  at  which  this  new  sunrise  line  crosses  the  great  circle 
path  (P^)  is  shown  in  Figure  II-3.   The  coordinates  of  this  point  can 
be  found  by  spherical  trigonometry: 

,      v  v  cos%(Ki    -   6<  o) 

tan^a   +  d)    =   tan?g(Ao   -    Ai) —, ~ =T" 

1  L     cos^C^j^   +^2) 

1   .u               1  ,                    n    sin%(a  +  d) 
cot%  f     =  tan%(0<1   -  **2> — 


sin%(a  -  d) 


sin&d   = 


sin  } 
sin  f 
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A'3      *A,  X 


A, 


Figure  II-3.      Revised   position  of   sunrise  point 
From  Figure   II -4: 

cos  p    =   sin  Xx   cos(  \3  -   \  x) 


cot 


^  cotAd  cos  ^3  -   sinAg  cos>C 

sin^O 

A3  =    A3  +   9 


sin  /3  \  =  sin  /$  3  cos^d   +  cos/S  3  sin^d   cosyO 


Pole 


<SO°-03 


Figure   II -4.      Sunrise   position  coordinates 
The  coordinates  of  the  point  of   sunset    (/3L>    A4)    can  De   found   in 
a  similar  manner  from  Figures   II-5   and   II-6: 

A4  =    A3   -    180° 
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A 

Ad' 
cot  e' 


=    "03 

=  -Ad 

cot  Ad*    cos^^  +  sin^^  coSiO 


sin 


sin/S^  -  sin^^  cos  Ad  -  cos/J^  sin  Ad  cos^ 


Figure  II-5.      Revised   position  of   sunset  point 


potr 

Figure  II-6.   Sunset  position  coordinates 
It  should  be  noted  here  that  if  it  is  assumed  that  (3 \   and  Ad' 
are  negative  quantities,  then  9'  will  also  be  negative  and  the  resulting 
longitude  correction  will  be  in  the  opposite  direction  from  the  sunrise 
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longitude  correction. 

Once  the  coordinates  of  the  four  required  points  are  determined 
the  amount  of  sunlight  can  easily  be  found  with  the  aid  of  Figure  II-7. 
This  shows  the  four  possible  situations  that  can  exist  on  the  great 
circle  path  of  propagation.   Points  (1)  and  (2)  are  the  transmitter  and 
receiver  positions,  (3)  is  the  sunrise  point,  and  (4)  is  the  sunset 
point.   The  amount  of  daylight  on  the  short  path  for  case  A  is  the  dis- 
tance from  (3)  to  (2);  the  amount  of  daylight  on  the  long  path  for  case 
A  is  the  distance  from  (4)  to  (2) ,  and  so  on. 

R 


case  A 


case  B 


CM  55 

case  C 


case  D 


Figure  II-7.   Possible  positions  of  sunrise  and  sunset 
on  great  circle  path 

The  actual  program  is  shown  on  pages  116  through  120.   The  geodesic 
program  was  converted  to  a  subroutine  for  convenience  here.   In  the  pro- 
gram the  declination  of  the  sun  as  a  function  of  the  time  of  the  year 
is  approximated  by  two  sine  waves  of  different  periods;  one  being  used 
from  the  1st  to  the  80th  and  the  267th  to  the  365th  days,  and  the  other 
being  used  from  the  81st  to  the  266th  day.   Data  from  the  Nautical  Al- 
manac Ql^l  was  used  as  a  guide  in  setting  up  this  approximation. 

The  program  also  includes  a  correction  factor  for  the  Greenwich 
Hour  Angle,  which  was  also  found  necessary  from  data  in  the  Nautical 
Almanac.   This  correction  factor  has  a  maximum  value  of  4.1  degrees 
(this  changes  ^\~  i-n  Figure  II-l)  and  is  again  approximated  by  sine 
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waves  of  various  periods  and  amplitudes  (see  statements  311,  511,  711, 
and  812  in  the  program) . 

Geocentric  latitudes,  obtained  from  the  geodesic  program,  are  used 
so  that  no  error  is  involved  due  to  the  ellipticity  of  the  earth.  The 
radius  of  the  earth  is  not  even  used  in  this  program,  since  the  distances 
between  points  are  expressed  as  arc  distances  and  the  amount  of  daylight 
is  expressed  as  a  percentage  of  the  long  or  short  path. 

The  only  inputs  required  are  the  time  of  day  and  the  day  of  the 
year.  The  latitude  and  longitude  of  the  two  points  of  interest  and 
the  description  of  the  great  circle  path  through  these  points  are  ob- 
tained from  the  geodesic  subroutine. 

In  the  paper  by  Brady  and  CrombieC223  three  basic  simplifications 
were  made  which  could  result  in  some  error.  These  were  the  assumptions 
that: 

1)  the  sun  is  a  point  source, 

2)  atmospheric  refraction  can  be  neglected, 

3)  the  earth  is  spherical. 

The  third  assumption  does  not  introduce  any  error  in  the  program  devel- 
oped here  because  geocentric  latitudes  were  used.   The  other  two  could, 
according  to  Brady  and  Crombie,  result  in  errors  of  about  16  minutes 
and  34  minutes  of  arc  respectively  in  the  quantity  £d  (Figure  II-3). 
These  errors,  although  small,  could  cause  considerable  error  on  high 
latitude  great  circle  paths. 

This  program  is  not  valid  for  the  case  in  which  the  highest  latitude 
on  the  path  is  greater  than  the  quantity  (90°  -  tXg) .   This  causes  a  por- 
tion of  the  path  to  be  continually  in  daylight  or  darkness   (e.g.,  polar 
paths  during  the  winter  or  summer)  and  the  method  used  here  breaks  down. 
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APPENDIX  III 
LOOP  ANTENNA  CALCULATIONS 
Consider  Figure  III-l  below.   This  is  a  schematic  representation 
of  a  loop  antenna  showing  a  voltage  generator  in  each  vertical  side  of 
the  loop  to  indicate  the  voltage  induced  by  the  passing  wave.   The  z- 
axis  is  taken  as  the  vertical  axis  and  the  xy  plane  the  horizontal  plane 
with  respect  to  the  earth.   The  voltage  ER  appearing  across  the  tuning 


Figure  III-l.  Loop  Antenna 
capacitor  C  is  the  signal  produced  by  the  antenna  and  is  found  by  phasor 
addition  to  be  ER  =  E2  -  E-^.   The  symbol  ~  here  denotes  a  phasor  quan- 
tity, while  a  letter  without  this  indicates  magnitude  only. 

In  order  to  find  the  effect  on  the  antenna  output  of  various  di- 
rections of  arrival  of  the  signal,  assume  the  antenna  to  be  rotatable 
about  the  z-axis.   The  angle  <X  is  defined  as  the  angle  between  the 
normal  to  the  loop  and  the  direction  of  arrival  of  the  signal,  as  shown 
in  Figure  III -2. 


DlJ2E.CnoM  OF 


ARRIVAL 


Figure  III-2.   Loop  Geometry 
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The  angle   varies   from  0°   to    360°   and  taken  to    increase  positively  as   the 
loop  rotates   clockwise.      The  bar   shown  on  the  loop   in  Figure  III-2   is 
for  reference  and  corresponds  to   side   2  of  the   loop  as   seen   in  Figure 
III-l. 

The  phasor  diagram  of  the  loop  voltages  is  given  in  Figure  III-3. 
The  angle  if  represents  the  phase  of  the  antenna  output  signal,  where 
(p  =  0  when  the  loop  is  oriented  broadside  to  the  direction  of  arrival 
of  the  wavefront. 


~ff 


Figure  III -3.   Phasor  Diagram 
The  angle  9  in  Figure  III- 3  represents  the  difference  in  phase  be- 
tween the  voltage  induced  in  side  2  and  that  in  side  1  of  the  loop.   It 
is  determined  by  the  distance  d  between  the  two  vertical  sides,  and  is 
given  by 


0  = 


2TTd 


sinOC  radians,  where 


A  =  free-space  wavelength  of  the  received  signal 
(same  units  as  d). 

Since  E-,  =  E2,  the  phasor  diagram  forms  an  isosceles  triangle  and  there- 
fore lf>  =  %9.  The  phase  angle  ip   is  then 


<P   = 


-     TT  d 
X 


sin  0C     radians, 
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The  resultant  ER  is  also  found  readily  from  the  phasor  diagram. 

Using  the  Law  of  Cosines,  one  has 

2     ?  2 

ER  =  E^  +  E2  -  2E1E2  cos  9 

=  2E^  -  2eJ  cos  e 

=  2E2  (1  -  cos  9) 

ER  =  E1  *[Z   (1  -  cos  6) 

=  Ei  V2  (Vfsin  — ) 

=  2EX  sin  <p  . 

It  is  convenient  to  normalize  by  setting  Ei  =  1.   One  can  then  express 
the  magnitude  of  the  relative  output  E   .  in  db  below  ER   : 


Eout  =  20  loSlO  "E 


% 


Rmax 


-   2°  log10  (  2  sin  I     >  '  Mhere  Yo  "  ^max 


o^  -1      f   sin  <0  N 
=  2°  loSlO  (  sin  fQ  )  ' 


123 


APPENDIX  IV 
CARDIOID   ANTENNA  CALCULATIONS 
As    in  Appendix  III,    assume  the  direction  of  travel   of  the   incoming 
signal    fixed   and  the   loop  antenna  rotatable,   with  the  angle   OC   defined 
as   before;    see  Figure   IV- 1.      In  addition,   the  sense  antenna   is  placed 
at   some  arbitrary  position,    forming  the   angle     OC'    as   shown.      The  distance 
d     of   the  sense   antenna  from  the  loop  midpoint   is   non-critical,    so   long 
as    it    is   short  with  respect   to   a  wavelength. 


SEMSE. 

AMTE.V1WA. 


DIRECTION   OF 
SIG,UAL  TRAVEL 


Figure   IV- 1.      Card io id  Geometry. 

Because  of  the  180°  phase  reversal  of   the   loop  antenna  when 
OC    =  0°   and  180°,    it   is   convenient  to  analyze  the  cardioid  system   in 
two  parts,   with  0°  <  OC  <   180°   and   with   180°<  OC  <    360°    (or  equivalently, 
-180°<  OC  <  0°)  .      As   in   Appendix   III,   the   loop  phasor  resultant   ER  is 
given  by  E2   -   En    and  phase    \$    is  8/2.      This    is   illustrated   in  the   phasor 
diagram  of  Figure   IV-2.      In  this   diagram,    both  6  and    if   are  taken  to  be 
negative  angles. 

For  proper  adjustment  of   the  cardioid   pattern,   the  magnitude  and 
phase  of  the  signal   from  the   sense  antenna  must  be  made  equal,    respec- 
tively,  to  the  magnitude  and   phase  of  the  signal   from  the  loop  antenna 
when  the   loop   is   positioned   for  maximum  response    (  0L  =   90°).      This 
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Figure  IV-2.  Loop  Phasor  Diagram 

condition  is  shown  in  Figure  IV-3,  where  Eg  is  the  magnitude  (after  the 

necessary  amplification  and  attenuation)  of  the  sense  antenna  signal. 

The  angle  ip Q   is  the  largest  (negative)  value  of  i^  .   Once  set,  the 

magnitude  and  phase  of  E  ,  of  course,  remain  constant  as  the  loop  is 

s 

rotated. 


Figure  IV-3.   Cardioid  Adjustment 
The  phasor  diagram  for  some  value  of  OL  between  0°  and  180°,  other 
than  <X  =  90°,  is  shown  below  in  Figure  IV-4.   The  overall  signal  out- 
put of  the  cardioid  unit  is  given  by  E0/Y°  (note  that  T  is  negative  and 
its  magnitude  less  than  that  of  ^0)«   This  is  obtained  by  the  phasor 
relationship  E  =  ER  +  Eg,  accomplished  in  the  cardioid  unit  which  mixes 
the  loop  and  sense  antenna  signals. 
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Figure  IV-4.   Card io id  Phasor  Diagram 
As  in  Appendix  III,  the  angles  (fi ,  <p  ,  0,  and  9Q  may  now  be  ex- 
pressed in  terms  of  the  angle  of  loop  orientation  01  ,  as  follows: 


9  =  ~d  sin  a   x  360°    (degrees) 


-d  sin  OL 


x  180O 


e„  =  9 


<Pn   =  <P 


-d 
-d 


max 


A 


x  360° 


x  180° 


The  objective  is  to  find  f  and  E  in  terms  of  OL  .   At  this  point, 
it  is  convenient  to  derive  an  expression  for  Er  as  a  function  of  E  . 
This  is  done  as  follows,  referring  to  Figure  IV-2  and  using  the  Law  of 
Cosines: 

E|  =  E^  +  E^  -  2E1E2  cos  6 


=   2E|    (1   -   cos  9),        since   Ex   =   E2. 
/,    ER  =    El  VT  s  in  9 


9   =  sin  9,        since  9   is   small 


Er 


=   9  VT. 
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E 


Rmax 


=  e0VT  = 


er    _ 


Figure   IV-5   is   an  expanded  view  of   the  phasor  diagram  of   Figure 
IV-4,   drawn  on  an  exaggerated  scale  to  make  angular  relationships 
clearer.      In  this  diagram  and  the  following  equations,   the   fact  that 


Figure   IV-5.      Angular   Relationships 
$    and   iQ      are  negative   is  taken  into  account.      In  the  triangle  of  Fig- 
ure  IV-5,  the  two   sides   ER  and  Eg   and  the  angle  between  them  are  known, 
so   that    X    may  be  found  by  geometry.      One  can  then  determine      i   =  if  -  ^ 
which  is   the  phase  angle  of  E   .      Solving  the  triangle, 


^    =   90°   - 


18QQ   +   y0-   if 


-1 


Es    -   ER 
Es   +    % 


L      s 


eot(   180°  t   »n  -  <f    )1 


=   90°   -   90°   +     S&gJEa      *  tan"1 


■  1  - 

9 
©o 

_  1  + 

e 

l(J    -   (P0 
cot  (90° 7— I ) 
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■'■  s  ■ 


1   -   1°     ♦  tan"1 


1  -  x  lf>    -  (*>„ 

tan 


1   +  x 


where  x  = 


Since    ({?  and    "C    can  be   found  for  any  given  value  of   01  ,   T    is  uniquely 
determined.      It    is   a  negative  angle   and    (f0<'>r<  0°   for  0°  <CX<180°. 

Determining  the  relative  amplitude  response   for  these  values  of    0L 
is  very  straightforward.      Refer  to   Figure   IV-5  and   use  the  Law  of   Cosines 
as   follows: 

Eg   =   E2    +  E'2    -    2EREg    cos    [l80°    +    (  lpQ   -  (p )] 

=   E2   +   E2    +   2ERES    cos    (<p      -  (J?) 


_   ^2 


R        "s 

E 


=    E 


R 


(-^-)2   +   1   +   2  -J£-     cos    (vpQ   -  ^  ) 

s  s 


e 


Eo  =  EsVl  +   C-|— ^"+  2"^—  cos^o  "  *P  ) 


=  V*   +  x2   +   2xcos    ( ip  Q   -  (f  ) 


One  may  now  normalize  by   setting   Eg    =   1   and  taking  201og-iQEo   to  obtain 
relative  output    in  db,    with  zero  db   occurring  at   0L    =0°   and    180°: 


lQ   =   20   log1QVl   +  x2   +   2x  cos(lpQ   -  Cf  ) 


db. 


The  procedures   for  obtaining  cardioid  phase  and  amplitude   for  loop 
positions    180°  <  QL  <  360°   are  quite   similar  to   the  methods   given  above. 
Figures   IV- 6   and   IV-7   below  are  self-explanatory  and  are   analogous   to 
Figures   IV-2   and  IV-4,   respectively,   discussed   above. 

Angles    are  similar   to  those  determined  before,   the  principal   dif- 
ference being  that    180°   is   now  added  to    tf    due  to  the  difference   in   loop 
position.      Note  that   all  of   the  following  angles   are  now  positive,   while 
"Y    remains   negative: 
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Figure   IV-6.      Loop  Phasors    (180°  <  0L<  360°) 


Figure   IV-7.      Cardioid  Phasors    (180°  <  a  <  360°) 


9  =   -  _!_  sin  a     x  360° 

=  180° —  sinQC   x  180° 


%   =  -|-  x  180° 


90  =  -~  x  360° 


As  before,    "C    is   found  by  trigonometry  as   follows: 

(<Po   +    ^    ~   180°)  -1 

XJ     =  90° 2-2 I +  tan   x 


Es   -  ER 
Es   +  ER 


cot    (■ 


Vo+I^-   180° 


=  180o  .   *£_!*_  +  tan-l 


x-1   M         4o   +   $ 

— t  tan 

x+1  2 


129 


IT  I  =  180°  -   [(  q>  -  180°)  +  \] 


.  T   =  _ 


180°  -  [(  If   -  180°)  +  \][ 


=  -180°  +  If  -  180°  +  TJ 

T  =  IP  +  ^  -  360°. 

Relative  amplitude  is  found  as  before  using  the  Law  of  Cosines,  and 
referring  to  Figure  IV-7: 

Eo  =  Es  +  ER  "  2ERES  cos  [180°  -  (^o  +  (f>] 
.'.  -jp-  =  Vl  +  x2  +  2x  cos  (lfQ  +  tf)  . 
Normalizing  and  expressing  the  result  in  db,  one  obtains: 

EQ  =  20  log10  V  1  +  x2  +  2x  cos(^Q  +  <^  )   db. 
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